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I. On the Mechaiiical Conditions of a Swarm of Meteorites, and on Theories of 

Cosmogony, 

lly G. H. Dabwik, LL.D.^ F.E.S.^ Felloiv of Trinity College and Plumian Professor 

in the University of Cambridge, 

Eeceiyed July 12, — Eead Noveinber 15, 1888. 

Mr. LocKiYER writes in his interesting paper on Meteorites'^ as follows : — 

'' The brighter lines in spiral nebulae, and in those in which a rotation has been set 
up, are in all probability due to streams of meteorites with irregular motions out of 
the main streams, in which the coUisions would be almost nil. It lias already been 
suggested by Professor G. Daiiwin {' Nature,' vol. 31, 1884-5, p. 25)— using the gaseous 
hypothesis — that in such nebulae 'the great mass of the gas is non-luminous; the 
luminosity being au evidence of condensation along lines of low velocity, according to 
a well known hydrodynamical law. From this point of view, the visible nebula may 
be regarded as a luminous diagram of its own stream-lines. ' '' 

The whole of Mr. Lockyer's paper, and especially this passage in it, leads me to 
make a suggestion for the reconciliation of two apparently divergent theories of the 
origin of planetary systems. 

The nebular hypothesis depends essentially on the idea that the primitive nebula is 
a rotating mass of fluid, which at successive epochs becomes unstable from excess ol 
rotation, and sheds a riug from the equatorial region. 

The researches of EocHEt (apparently but little known in this country) have 
imparted to this theory a precision which was wanting in Laplace's original 
exposition, and have rendered the explanation of the origin of the planets more 
perfect. 

But notwithstanding the high probability that some theory of the kind is true,J the 
acceptance of the nebular hypothesis presents great difficulties. 

Sir William Thomson long ago expressed to roe his opinion that the most pro- 
bable origin of the planets was through a gradual accretion of meteoric matter, and 

^ ' Nature,' Nov. 17, 1887. The paper itself is in ' Roj. Soc. Proc.,' Nov. 15, 1887 (No. 259, p. 117). 
t ' Montpellier, Acad. Sci. Mem.' 

t [The very remarkable photograph of the nebula in Andromeda, exhibited to the Royal Astronomical 
Society by Mr. Isaac Boberts on December 6, 1888, affords something like a proof of the substantial 
truth of the nebular hjpoihesis.-~-a. H. D. December 19, 1888.] 
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the reseai-ches of Mr. Lockyek aflford actual evidence in favour of the abundauce of 
meteorites in space. 

But the very essence of the nebular hypothesis is the conception of fluid pressure, 
since without* it the idea of a figure of equilibrium becomes inapplicable. Now, at 
first sight, the meteoric condition of matter seems absolutely inconsistent with a fluid 
pressure exercised by one part of the system on another. We thus seem driven either 
to the absolute rejection of the nebular hypothesis, or to deny that the meteoric con- 
dition was the immediate antecedent of the Sun and Planets. M. Faye has taken the 
former course, and accepts as a necessary consequence the formulation of a succession 
of events quite dififerent from that of the nebular hypothesis. ""'^ I cannot myself find 
that his theory is an improvement on that of Laplace, except in regard to the adop- 
tion of meteorites, for he has lost the conception of the figure of equilibrium of a 
rotating mass of fluid. 

The object of this paper is to point out that by a certain interpretation of the 
meteoric theory we may obtain a reconciliation of these two orders of ideas, and may 
hold that the origin of stellar and planetary systems is meteoric, whilst retaining the 
conception of fluid pressure. 

According to the kinetic theory of gases, fluid pressure is the average result of the 
impacts of molecules. If we imagine the molecules magnified until of the size of 
meteorites, their impacts will still, on a coarser scale, give a quasi-fluid pressure. I 
suggest then that the fluid pressure essential to the nebular hypothesis is, in fact, the 
resultant of countless impacts of meteorites. 

The problems of hydrodynamics could hardly be attacked with success, if we w^ere 
forced to vstart from the beginning and to consider the cannonade of molecules. But 
when once satisfied that the kinetic theory will give us a gas, which, in a space 
containing some millions of molecules, obeys all the laws of an ideal non-molecular 
gas filling all space, we may put the molecules out of sight and treat the gas as a 
plenum. 

In the same way, the dlflJculty of tracing the impacts of meteorites in detail is 
insuperable ; but, if we can find that such impacts give rise to a quasi-fluid pressure on 
a large scale, we may be able to trace out many results by treating an ideal plenum. 
Laplace's hypothesis implies such a plenum, and it is here maintained that this 
plenum is merely the idealisation of the impacts of meteorites. 

As a bare suggestion this view is worth but little, for its acceptance or rejection 
must turn entirely on numerical values, which can only be obtained by the considera- 
tion of some actual system. It is obvious that the solar system is the only one about 
which we have sufficient knowledge to afford a basis for discussion. This paper is 
accordingly devoted to a consideration of the mechanics of a swarm of meteorites, with 
special numerical application to the solar system. 

The investigation has entailed a considerable amount of mathematical analysis ; 

* ^Snr rOrigine du Monde/ Paris, Gauthibr-Yillars, 1884; * Annuaire pour 1' an 1885, Bureau des 
Longitudes,' p. 757. 
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there is, however, no analysis in §§ 1 and 2. The reader who only wishes to know 
the arguments and results, without a consideration of the mathematical details, is 
therefore recommended, after reading §§ 1 and 2, to pass on to the Summary. 

§ 1. On the Effective Elasticity of Meteorites in Collision, 

When two meteoric stones meet with planetary velocity, the stress between them 
during impact must generally be such that the limits of true elasticity are exceeded ; 
and it may be urged that a kinetic theory is inapplicable unless the colliding particles 
are highly elastic. It may, however, I think, be shown that the very greatness of 
the velocities will impart what virtually amounts to an elasticity of a high order of 
perfection. 

It appears, a priori, probable that, when two meteorites clash, a portion of the 
solid matter of each is volatilised, and Mr. Lockyer considers the spectroscopic 
evidence conclusive that it is so. There is, no doubt, enough energy liberated on 
impact to volatilise the w^hole of both bodies, but only a small portion of each stone 
will undergo this change. 

A rough numerical example will show the kind of quantities with which we are 
here dealing. 

It will appear hereafter that the mean velocity of a meteorite may be at the least 
about 5 kilometres a second; and, accordingly, the mean relative velocity of a pair 
would then be about 7 kilometres a second.^ Hence, if two stones, weighing a 
kilogramme, move each with a velocity of 3^ kilometres per second directly towards 
one another, the energy liberated at the moment of impact is 2 X ^ X 10^ (3| X 10^)^ 
or 12 X 10^^ ergs. 

Now Joule's equivalent is 4*2 X 10^ ergs; hence, the energy liberated is about 

3 million calories. 

It is quite uncertain how much of each stone would be volatilised ; but, if it were 
3 grammes, there w^ould be a million calories of energy applied to each gramme. 

The melting temperature of iron is about 1500 degrees Centigrade, and the mean 
specific heat of iron may be about y.t Hence, about 300 calories are required to raise 
a gramme of iron from absolute zero to melting point. I do not know the latent heat 
of the melting of iron, but for platinum it is 27, and the latent heat of volatilisation 
of mercury is 62. Hence, about 400 or 500 calories suffice to raise a gramme of iron 
from absolute zero to volatilisation. It is clear, then, that there is energy enough, 
not only to volatilise the iron, but also to render the gas incandescent ; and the same 
would be true if the mass operated on by the energy were 30 grammes instead of 3. 

It must necessarily be obscure as to how a small mass of solid matter can take up 
a very large amount of energy in a small fraction of a second, but spectroscopic 
evidence seems to show that it does so ; and, if so, we have what is virtually a violent 
explosive introduced between the two stones. 

^ If ^ be the velocity of mean square, t? \/ 2 is the square root of the mean square of relative velocity, 
f 'Physikalisch-Chemische Tabellen.' Landolt and Boknstein. 
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In. a direct collision each stone is probably shattered into fragments, like the 
splashes of lead when a bnllet hits an iron target. But direct collision must be 
a comparatively rare event. In glancing collisions the velocity of neither body is 
wholly arrested, the concentration of energy is not so enormous (although probably 
still sufficient to effect volatilisation), and, since the stones rub past one another, more 
time is allowed for the matter round the point of contact to take up the energy ; thus, 
the whole process of collision is much more intelligible. The nearest terrestrial 
analogy is when a cannon-ball bounds off the sea. In glancing collisions fracture will 
probably not be very frequent. 

From these arguments, it is probable that, when two meteorites meet, they attain 
an effective elasticity of a high order of perfection ; but there is, of course, some loss 
of energy at each collision. [It must, however, be admitted that on collision the 
deflection of path is rarely through a very large angle. But a succession of glancing 
collisions would be capable of reversing the path ; and, thus, the kinetic theory of 
meteorites may be taken as not differing materially from that of gases. "^ 

Perhaps the most serious difficulty in the whole theory arises from the fractures 
which must often occur. If they happen with great frequency, it would seem as if 
the whole swarm of meteorites would degrade into dust. We know, however, that 
meteorites of considerable size fall upon the Earth ; and, unless Mr. Lookyer has 
misinterpreted the spectroscopic evidence, the nebulae do now consist of meteorites. 
Hence, it would seem as if fracture was not of very frequent occurrence. It is easy 
to see that, if two bodies meet wdth a given velocity, the chance of fracture is much 
greater if they are large, and it is possible that the process of breaking up will go on 
only until a certain size, dependent on the velocity of agitation, is reached, and will 
then become comparatively unimportant. 

When the volatilised gases cool, they will condense into a metallic rain, and this 
may fuse with old meteorites whose surfaces are molten. A meteorite in that 
condition will certainly also pick up dust. Thus, there are processes in action 
tending to counteract subdivision by fracture and volatilisation. The mean size of 
meteorites probably depends on the balance between these opposite tendencies. If 
this is so, there will be some fractures, and some fusions, but the mean mass will 
change very slowly with the mean kinetic energy of agitation. This view is, at any 
rate, adopted in the paper as a working hypothesis. It was not, however, possible 
to take account of fracture and fusion in the mathematical investisfation, but the 
meteorites are treated as being of invariable mass. 

§ 2. On the Velocity of Agitation of Meteorites^ and on its Secular Change, 

The velocity with which the meteorites move is derived from their fall from a great 
distance towards a centre of aggregation. In other words^ the potential energy of 
their mutual attraction when widely dispersed becomes converted, at least partially, 

* Added FoY. 16, 1888. 
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into kinetic energy. When the condensation of a swaxm is just beginning, the mass of 
the aggregation towards which the meteorites fall is small ; and, thus, the new bodies 
arrive at the aggregation with small velocity. Hence, initially, the kinetic energy is 
small, and the volume of the sphere within which hydrostatic ideas are (if anywhere) 
applicable is also small. 

As more and more meteorites fall in that volume is enlarged, and the velocity with 
which they reach the aggregation is increased. Finally, the supply of meteorites in 
that part of space begins to fail, and the imperfect elasticity of the colliding bodies 
brings about a gradual contraction of the swarm. 

I do not now attempt to trace the whole history of a swarm, but the object of the 
paper is to examine its mechanical condition at an epoch when the supply of meteorites 
from outside has stopped, and when the velocities of agitation and distribution of 
meteorites in space have arranged themselves into a sub-permanent condition, only 
affected by secular changes. This examination will enable us to understand, at least 
roughly, the secular change in the velocity and in the distribution of the meteorites 
as the swarm contracts, and will throw light on other questions. 

I 3. Formulm for Mean Square of Velocity, Mean Free Path, and Interval between 

Collisions. 

We have to investigate whether, when the solar system consisted of a swarm of 
meteorites, the velocities and encounters could have been such that the mechanics of 
the system can be treated as subject to the laws of hydrodynamics. The formulae 
which form the basis of this discussion will now be considered. 

For the sake of simplicity, the meteorites will, in the first instance, be treated as 
spheres of uniform size. 

The sum of the masses of the meteorites is equal to that of the Sun, for the planets 
only contribute a negligible mass. If M^ be the Sun s mass, and m that of a meteorite, 
their number is MJm. 

If, at each encounter between two meteorites, there were no loss of energy, the sum 
of the kinetic energies of all the meteorites would be equal to the potential energy 
lost in the concentration of the swarm from a condition of infinite dispersion, until it 
possessed its actual arrangement. In such a computation the rotational energy of the 
system is negligible. 

Suppose the Sun's mass to be concentrated from infinite dispersion until it is 
arranged in the form of a homogeneous sphere of radius a and density p. Then let 
the sphere be cut up into as many equal spaces as there are meteorites, and let the 
matter in each space be concentrated into a meteorite. When the number of 
meteorites is large, the potential energy lost in the first process is very great compared 
with that lost in the subsequent partial condensation into meteorites.^ Thus, the 
energy lost in the partial condensation is negligible. 

* It depends, in fact, on the square of the ratio of the diameter 2a to the linear dimersion of one of 
the eqnal spaces. 
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If fjb be the attractlonal constant, the lost energy of condensation is well known to 
be ffiMQ^/a. But on the hypothesis that there is no loss of energy at each encounter, 
this must be equal to the sum of the kinetic energies of all the meteorites. If, there- 
fore, v^ be the mean square of velocity of a meteorite, we must have ^ MqV^ = ^iiM^ja, 
so that v^ =. ^{JL MJa. 

But homogeneity of density and uniformity of kinetic energy of agitation are 
impossible ; for the meteor-swarm must be much condensed towards its centre, so that 
we have largely underestimated the lost potential energy of the system. Also, the 
velocity of agitation must decrease towards the outside, or else the swarm would 
extend to infinity. Besides this, the partial conversion of molar into molecular energy, 
which must take place on each encounter, has been neglected. 

We shall see below reason for believing that throughout a large central volume the 
mean square of velocity of agitation is nearly uniform, and that outside of this region 
it falls off. 

Suppose, then, that If is the mass and a the radius of that portion of the swarm 
in which the square of velocity of agitation is uniform ; let /Vq^ be that square of velocity, 
and let it be defined by reference to the poteutial of M at distance a, so that 



v^ 



3 — m 



fiM 







f'~- « 



where /3 is a coefficient for which a numerical value will be found below. 

The square of velocity of agitation outside of the radius a is to be denoted by #, 
and subsequent investigation will be necessary to evaluate v^ in terms of v^^. 

If we denote by a^ the Earth's distance from the Sun, and by u^ the Earth's 
velocity in its orbit, we have 



Whence, 



Uq^ — • jM • ^2) 



If in any distribution of meteorites %v is the sum of the masses of all the meteorites 
in unit volume, or the density of the swarm at any point, and if X be that distance 
which is called in the kinetic theory of gases '' the mean distance between neighbouring 
molecules," w^e have 

•~ o».»..»...iy:l 

w ■ ^ ^ 

Now, the mean density of that part of the sw^arm in which the kinetic energy of 
agitation is constant being p, we have 

anci 
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Suppose that s m '^ the radius of the sphere of action" of a meteorite^ so that when 
two of them approach so that the distance between their centres is s there is a 
coUision. 

Let L and T be the mean free path and mean interval between collisions. Then, 
since the mean velocity is v \/(8/3 tt), we have, according to the kinetic theory of 
gases,"^ 

X = -^. r=z^,/f7r (7) 



TT^V^ 



V 



Then, on substitution from (4), (5), and (6), we have 
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■ ^0 V8" 
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and we have 
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(8) 



. . (9) 



(10) 



We now proceed to calculate Uq, Zq, Tq, and also 2^0//^, using the centimetre-gramme- 
second system of units. 

The Sun s mass may be taken as 315,511 times that of the Earth, and the Earth as 
6*14 X 10^'^ grammes t; hence 

M^ = 10^3'^^7i8 — 1-9372 X 10^^ grammes. 
The attractional constant and the Earth's mean distance from the Sun are 






a 







1'487 X 10'^ cm. 



* Meyer, ' Kinetisclie Theorie der Oase.' 

t Here and elsewhere I generally use Eyerett's ^ Units and Physical Constants.' 



8 PROFESSOR G. H. DARWIN ON THE MECHANICAL CONDITIONS OF 



With these vahies 



Mg = 106-46323 __ 2,905,600 cm. 'per sec. "^ 
I = 106-6812S = 4,800,600 cm. 






100-25366= 179334 see. V - (H) 

106-79204 __ 6^195^000 



^ 



The dimensions of l^ and Tq are not those of length and time ; but, if meteorites 
of 1 gramme mass, with sphere of action 1 centimetre, and '^ velocity of mean square ^' 
of agitation equal to the Earth's velocity in its orbit, have density of distribution 
equal to one-third of the mean density of the sphere M, then ?q, Tq will be the mean 
free path and time, as stated in centimetres and seconds. We may thus regard Zq, t^ 
as a length and time, provided care be taken in the subsequent use of the symbols to 
adhere to the c.g.s. system of units. 

§ 4. On the Equilihrmm of a Gas at Uniform Temperature in Concentric Spherical 

Layers under its own Gravitation, 

It is assumed provisionally that the conditions are satisfied which permit us to 
regard the swarm of meteorites as a quasi-gas, subject to the laws of hydrostatics. 

The solution of this problem, then, becomes a necessary preliminary to the discus- 
sion of the kinetic theory of meteorites. The equilibrium of a gas under its own 
gravitation has been ably discussed by Professor Emitter in one of his series of papers 
on gaseous planets."'"'' The intrinsic interest of the problem renders an independent 
solution valuable. Suppose, then, that a mass M-^ of gas is enclosed in a spherical 
envelope of radius a^ and is in equilibrium in concentric spherical layers. Let v^, 
the mean square of the velocity of agitation of the gaseous molecules, be defined by 
reference to the potential of the mass M\ at the radius a^, so that 



a^ 



where ^^ is a numerical coefficient, and /x is the attractional constant. 

Let p and w be the pressure and density of the gas at radius r, and k the modulus 

of elasticity, so that 

p :=^ lew, 



k = i v;^ = i A 



2 



a^ 



* " Untersuclmngen iiber die Holie der Atmosphare und die Constitution gasformiger Wieltkorper," 
'Wiedemann's Annalen ' (ISTew Series), vol. 16, 1882, p. 166. A very elegant solution of part of my 
problem has also been given by Mr. Gr. W. Hill in tbe ' Annals of Mathematics,' vol. 4, No. 1, p. 19 
(February, 1888). Mr. Hill's paper only readied my hands after my own calculations had been 
completed, and I therefore adhere to my own less elegant method. Mr. Hill has obviously not seen 
M. Eitter's papers. 
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Then the equation for the hydrostatic equilibrium of the gas is 



— 7- + 4 TT u wr^ dr = 0, 

w dr ' Jq 



(12) 



It is obvious that — is equal to the whole mass enclosed inside radius r, and 

/jLW dr -^ 

this relation will hold however the equation be transformed, provided we do not 

multiply the equation by any factor. 

In consequence of the relation between p and w this may be written 



k 



r^ ~ log w + 



AlTfJU 



wr^ dr 







= 0, 



If pi be the mean density of the mass M-^, we have 

4 TT /X = 



pi<^ ^iWpi 



Hence, we may write the equation (12) in the form 



ha^ 



Now, let 



'r d 
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a^ dr ^ ' /3i^ 



r^ w r^ , 
~ idr 

Pi (^i _ 






and the equation becomes 



w 



Pi 



= iA'>A 






0. 



= 0. 



• • • 



By differentiation we obtain the equation 



• (13) 



Vvt/Ji t«/-i 



0. 



(14) 



It is obvious from (13) that ^ fi^^ M^ dyjdx^ is the mass enclosed inside radius ajx^, 
and therefore -J ^^ dy^jdx^ is equal to unity when x=l. 

A general analytical solution of (14) does not seem, to be attainable, and recourse 
must be had to numerical processes. Although this is an equation of the second 
degree, and its general solution must involve two arbitrary constants, we shall see (as 
pointed out by M. Eitter) that the general solution, as applicable to our problem, 
may be deduced from one single numerical solution. M. Eitter proceeds by a 
graphical method, which he has worked with surprising accuracy. I shall therefore 
adopt an analytical and numerical method, which, although laborious, is susceptible of 
greater accuracy. 
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Whatever be the arrangeroent of the gas, the density at the centre must have some 
value. I therefore start with a central density co, corresponding to the value rj of y^, 
j5o that 

"^ = i^,^e\ (15) 



CO 

Pi 



For the sake of brevity the suffixes 1 will be now omitted from the various symbols, 
to be reaffixed later when they are required. 

A.t the centre, where x is infinite, dy/dx, d^yjdx^, &c., are all zero, and we put 

Let ^ = e^lx^, and let us assume 

NoW; the differential equation (14) to be satisfied is 



But 



Q(yl '1- :;:::;: — zz: ^e . 

dx^ x" 



x^'pl= -2. 3A,i + 4 . 5A,e - 6 . 7A,e + • ■ • 

dx^ 



and by expanding e'' we obtain 

- ^e— - ^ + A,e -{A + i A') e + (^3 + AA + A A,^) t 
- (A, + A,A, + i ^3^ + 4 A,U, + ^hr A') ^ 
+ {A, + A,A, + A,A, + ^A,U, + ^AA' + 2-3^1^43+ -^rs^^-T ^i') i' 

By equating coefficients in these two series, I find 



■^1 — 6 ? -^2 — 12 0? ^-3 —1890? ^' 



61 



4 — 1,63 2,960 ? 

A 629 4 — 3407383 ... ^p 

^5 ■— 224,532,000 ' ^6 ~ 156x10^ ' ^^'' 

and 

log A^ = 9-2218487, log A^ = 7*9208188, log A^ = 6*7235382, 

log A^ = 5-5723543, log ^5 = 4*4473723, log Aq = 3-3392964 ; 

whence, by extrapolation, 

log ^7 = 2-243, log ^8 = 1*13. 

In M. RiTTEn's paper, already referred to, he takes a certain function u as equal to 
1-031, when the radius is unity. Now, Eitter's function u is equal in my notation to 
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^ {iv/p) -r- ^ fi^a^/r^ or J- e^/x^. It follows, therefore, that Ritter takes the surface 
value of 2/ = log^ 2*062. But he intends the central density to be 100 times the 
surface density; hence, to take the same solution, we must have e^ eqiial to 100 
times 2*062. Therefore, his value of 77 would be 

rj = log, 100 + log, 2*062, 

or 

77 = 5*3288465. 

Kow, as I want to make a comparison between my solution ^nd his, I start with 
this value of 77. The only object attained by the choice of this particular value i^ 
that the two solutions become easily comparable. It will be seen below that the 
value of 77 does not make the central density exactly 100 times the surface density, 
but only satisfies that condition approximately. In Hitter's graphical treatment of 
the problem this value 100 is the exact datum, whereas in my method we start with 
an exact value of 77, and proceed to find the ratio of central to surface density. 
. With this value of 77 (whence log^o e'' = 2*3142888) I find the following series 
for y : — 

^0000.^^ 34-3667 , 354-321 4638-79 , 67,532*0 1,044,280 
y = 5*3288465 — H ~ — + — ^-r — -^ 

t^ lA/ tO I/O iJU 

I 1^^789,00 ^ 2-77 x 10^ 4-4 x 10^ _ ^ 

and, by differentiation, the series for dy/dx is obvious. 

. This series will be very accurate from ir= 00 to about cc = 8. Thus, when r/a = '1, 

or iT = 10, we have 

2/= 5*016558, f^= + -0568910, 

ax 

and even the far less convergent series for d^y/dx^ gives — '0150891, agreeing with 
-- e^/x"^ to the last place of decimals. When r/a = '125, or x = 8, we have^ 

V= 4*863925, ^- = *101168, 
whence, 

f f = - -031624, 

with y correct to four, and probably to five, places of decimals, and dy/dx probably 
correct to four places of decimals. This is amply sufficient for our purpose. Indeed, 
accuracy of this order would be altogether pedantic, were it not that the errors 
accumulate. 

* Even when x = b,l find from tliivS series y = 4*342, wMch lies very near to y = 4-332, found below. 
But the series for dyjdQi is useless. 
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We cannot, then, rely on this method of procedure beyond the region included 
tetween a? = oo and £C = 8, and must now make a new departure. 
Since 



dw^ 






log ( — ■ 



d?y\ 
dx' 



y — 4: log X ; 



therefore, 



d^y d^y fdy 4t 

dm^ dx^ \dx x 



% ^ 9 • 9 '^ ^ 



Now, let 



A 



n 



1 d^'y 



■ (ir) 



where, after differentiation, x is put equal to c, a constant. 
Then (17) may be written — 



Now, it is clear that 



A. = pjA,-*' 



dfA„ n + pi , 



• • 



. . . . (18) 



Hence, differentiating (18) n — 3 times, we have 



nl . 



q = n--z ^-.31 ^In — q — ll 



= t 



q-O 



n — 3 — ql ql 3! 



2! 



JUL 



q -\- 11 



n 



.,_4'----A« + (-)' 



4(gl) 



C2 



+ 1 r> 



or 



A — 



■n 



n,n^l.n-d,^o ^^''^'^'^) ^'^^'i'-^'^) ^--^-^[^^ 



cr 



or 



A — 



n.n — 1 .n — 6 ' 



(»-2)^„_3+(-) 



n 



4 



n — '^ 



6' 



+ 3 . 2^; 






How, if, for a given value of a^, viz., c, we know y, or J[q, and dyjdx, or ^j, then we 
can compute Ac^ from the formula — -le^^C^^; and, by the formula (18), viz. : — 



A 



3 



o . /j . X 



2AA 



2 



^1 



4" 



A^ may be computed. 

Afterwards^ A^, A^, &c., may be computed by successive applications of (19). 

being so^ 



This 
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y = Aq + Ai{x — c) + A^{x — c)^ -^ . . , , . . . . (20) 
^ = A^^ 2A^ (a? — c) + 3^3 {x — cf + . . , . 

In these series x may have any value, provided the series converges adequately. 
The convergence may be much improved by an artifice, which, however, I unfortu- 
nately did not discover until most of the computations were completed. Let us add 
and subtract log 2x^ on the right-hand side of (20). 

Now, 



logv 2x^ = log^ 2c^ + 2 log, 



1 + 



u/ ~~" C 



— lOg^ ZC -f Z — -2 -t- 3 



If, then, we write 

2 2 2 

^0 = ^0 — log, 2C^ 5i = ^i — ^, B^ = A^ + j^^r 53 = ^3— —&c., 

we have 

2/=log.2a?^ + Bo + -Si(^-^)+52(^-c)'+ • •*. • • • (21) 

a more convergent series than that with the ^'s. 

The simplest way of computing the B's appears to be by first computing the A's, 
The process for obtaining the numerical solution is then as follows : — 
We have the values of y, dy/dx, ^dy^/dx^ when cc = 8, that is to say, of ^q, A^, A^ 
w^hen c = 8. From these the successive ^'s and jB's are computed, and the resulting 
series gives the values of y and dy/dx when i^ is 5 or r = '2. Starting from this 
point a new series gives the result when r = '3, another series gives the values for 
r = -4, and so on. Later in the calculation several values may be computed from one 
formula.^ 

When the computation has been carried out to r = a, we have reached the 
spherical envelope, but that envelope may be replaced by another at any more remote 
distance from the centre. Thus, the integration may be pursued for values of x less 
than unity, and when the lower limit is zero the envelope is at infinity. 
If we write 

log u r= Bg + Bi{x — c) -{^ B^{x '■— cY "{-,,. , 



* If the series be carried as far as Bg, several steps may be included in one series. For example, the 
first series, when c = S, may be pushed even as far as r = '4 without serious error; for it gives j=: 2*960 
instead of the true value 2-965^ and dyjdx -= -957, instead of the trno value '944. I have not, however, 
been satisfied with this degree of accuracj^. 
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we have 

and 

w 2 /3^a^ 



9 Q 



,3 



• u. 



But it may easily be seen that 2y8%^/9r^ is a particular solution of the problem ; 
hence, u is a factor by which the particular solution is to be multiplied to obtain the 
general solution. The function u is given by 

^ ~ "^ ^ ~ ""* ^^" ~(M' ........ (22) 



A table of the values of u is given below, showing how the general solution shades 
off into the particular solution. This function, k., is also tabulated by Ritter, and I 
made use of its value, when cc = 1, to determine the value of tj, with which the 
integration is to begin. I find, however (see Table I.), that, when x = 1^ u= 1*0063,- 
in place of 1'031, as given by him. 

The last row in Table I. gives the ratio of the central density co to w, the density 
at the distance r ; this ratio is equal to e'^'^K 
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It will be noticed that u rises from zero to a maximum of about 1*66, falls to a 
minimum of about '82, and then rises to unity. 

Since -J jS^^ dyjdx^ = 1 when x^ = 1, we have ^ ^^^ = 1/2-4087 = -4152. 

M. EiTTEE has '4143 for this constant, which he calls m. 

It appears from the Table that the density at the centre is 102^ times as great as 
that where r = %. M. Hitter's solution is intended to make that ratio exactly 100, 
but this solution shows that we ought to have started with a slightly different value 
of 7) to obtain that result. 

In the general solution of the differential equation d^y/dx^ = — e^/x"^ the two 
arbitrary constants may be taken to be the values of y and dyjdx when x is infinite. 
Now, we have taken arbitrarily y = 5'329 when x is infinite, and the physical 
conditions of the problem imply that dy/dx is zero when x is infinite. For if dy/dx 
had any positive or negative value different from zero, it would mean that at the 
centre there was a " nucleus of infinitely small dimensions, but of finite positive or 
negative mass. Now, a^ is that distance from the centre at which the density is 
1/102*45 of the central density ; hence, we may regard a^ as the arbitrary constant of 
the solution. Whatever be the elasticity of the gas, we may always take as our unit 
of length that distance from the centre of the nebula at which the density has fallen 
to 1/102*45 of its central value. Hence, the above table gives the general solution of 
the problem, subject, however, to the condition that there is no central nucleus. 

If we view the nebula from a very great distance, a^ appears very small, and thus 
the solution of the problem becomes y = log 2x^. It is easy to verify that this is a 
particular algebraic solution of the differential equation, as is pointed out by Ritter 
in his paper. ^ I found this solution very useful in a preliminary consideration of the 
problem treated in this paper. 

The next point which we have to consider is the form which the solution will take, 
if, instead of taking a^ as the unit of length, we take any other value. 

The density at any distance and the elasticity are to remain unchanged, but are to 
be referred to new constants. 

Thus, w, r, v^ remain unchanged, but are to be referred to M, p^ /3^, a, instead of to 

Now, since w remains unchanged, 

and, since v^ remains unchanged, 
Also 

a a 

* I have made use of this solution in a paper in tlie ^ ProceedingB of tlie Eoyal Society,' Dec. S, 1883, 
and it has also been referred to in a paper by Sir W. Thomson, 'Phil. Mag./ vol. 23, p. 287. Sir W. 
Thomson's paper covers mnch the same gronnd as some of M.^Ritter's earlier papers, bnt was written 
by him independently and in ignorance of them. 
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From these relations it is clear that 



2/ = 2/1 - 2 log ^ , 
and 

Then, since ^ ^^ dyjdx = 1, when x= 1^ and since dy = dy^ and dx = dxi aja-^^ it 
follows that 



^2 -^ _____ — --— ^ when Xi=^-^» ...... (23) 



1 02 



JU-\ CviJ-i I CLOut (jj 



This relationship has been already used for determining ^^. 
It is obvious also that 

p x^ dy^ldx^, when x^^a^jct 

pi x-^ dy^ldxi, when Xi == 1 
Therefore, 



W ^6^1 



/> sc^dyijdxij wlien ^^^^ = ^^/(^ 



• (I 



• • (24) 



If Wq be the density when r =^ a, we have 



^^0 h^\ when x.=:aJ a 1 x^d^yjdx^^ , - . , . 

— = - — ' — Q-r-T^ ^^ ==""+• -~r^T~^ J when x^ = ai/a. . (25) 

p x{ dy^dx^ "^ dy^dx-^ ^ ^' ^ ^ 



If Pq be the pressure when r = a, we have 



Po = i^'^^0 = W(^y ■ ^ J 



If, therefore, we write P = ^irjxa^p^^ 

-^0 _ 1 /D2 ^0 d y^dx-^ , , /r»r»\ 

P - 3/3^ y - - (^^^^ , when a;^ = a^/a (26) 

By (26) we are able to find how the pressure on an envelope of given radius a 
varies with the variation of the temperature of a given mass M of gas contained in it. 
By means of the formulae (23), (25), (26), we are now able to obtain from the original 
solution any number of other ones ; for, after the changes have been effected in the 
notation, we may proceed to magnify or diminish all the various values of a until they 
are of one size, and we shall thus obtain the solution for a gas at any temperature 
whatever. 

I shall now proceed to give a table of results when the standard radius a, which 
may be conveniently called the boundary, is placed successively infinitely near the 
centre, where r = X a^, at r = '1 X a^, r = '2 X a^, and so on. The first line of 
entries gives the various values of ^yS^ (computed from (23)), on which the elasticity 
of the gas depends; the second line gives ivjp (computed from (25)), or the ratio of 

MDCCCLXXXTX. — A. D 
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the boundary density to the mean density of all inside of it ; the third line gives 

Po/P (computed from. (26)), by which we trace the variations of pressure at the 

boundary. 

Table II. 



Value of a by ref ex^ence 1 ^ 
to former solution J a^ 





•1 


-2 


-3 


•4 


•5 


•6 


•6264 


•7 


•8 


i^^ 


-iP' 


r -^ 1 




GO 


1-7577 


-6741 


-4826 


-4236 


•4031 


•3974 


•3972 


•3984 


•4027 


juMja 


Xi dyildxi_ 


P 




— --x-^ d^y^jdx^ 




1-0000 


•8841 
4-662 


-6838 


-5333 


•4383 


•3791 


•3410 


1 

3 


•3158 


•2986 


dyjdxi 


F 


— d^y-yldxi^~ 
[dyildxiY _ 




CO 


1-383 


•772 


•557 


•458 


•407 


•397 


•377 


•361 



Yalue of a loj reference 1 '^ _ 
to former solution j a^ 


•9 


1-0 


1^25 


1-5 


2-0 


2-5 


3-0 


00 


3^ 


= i/3^ 


- ^ 1 




•4085 


•4152 


-4325 


•449 


•476 


•487 


•497 


1 

2 


jLiMla 


x-^ dy^jdxY_ 


Wq 

P 


— -|a?2 d^yjjdxi^' 
dy^jdx^ 




•2867 


•2785 


•2676 


•264 


-267 


•269 


-273 


1 
3 




P(^ 
F 


"__ — dhj^jdx^^ 
[dyjdx^j'^ _ 


— 


•351 


•347 


•347 


-356 


-382 


•392 


•406 


1 

2 



The minimum value of w^j p occurs when aja^ =1*6 very nearly, for, when 
a/ai=l-4, 1-5, 1-6, I find wjp= '2&521, '26437, '26425 respectively.^ When 
TJa^ = 1*6, y^=z ~ -38435 and dyjdx^ == 3*5180. The minimum value of pJP 
occurs when aja^ = I'l very nearly, for, when a/% = I'O, 11, 1*2, I find p^jP = '3469 
'3455, '3462 respectively. 

When wjp is a minimum, the density at the centre is 381 times that at the 
boundary, and, when pJP is a minimum, the density at the centre is 129 times that at 
the boundary. M. Hitter finds the pressure to be a minimum when this ratio is 258, 
instead of 129. As this corresponds to aja^^^ 1*5, this discrepancy between our 
solutions is not so large as might be expected from the great discrepancy between 
these results, and I cannot but think that my result is more accurate than his. 

The minimum value of -g-/3^ occurs when a/a^ =: '6264, and its value is '39723. 
This value makes the surface density exactly one-third of the mean density, for -g- ^^ is a 
minimum when x^ dy^jclx^ is a maximum, and this occurs when x^d^yjdx-^ -f- dy^ldx^ ^^ ; 
and, when tliis relationship is satisfied, lo^/p = ^. 

It is interesting to note that in this case /S'^ is very nearly equal to |, so that the 



* Mr. Hill finds that fclie mioimnm value of w^/p approximates to -pg-, or -2667. The agreement 
between our results is satisfactory. 
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total internal kinetic energy of agitation of the sphere of gas at minimum temperature 
limited by the radius a is ^ (f /xlf/a) ilf = ^ [juM^a very nearly. Now/ the energy 
lost in the concentration of a homogeneous sphere M from a condition of infinite 
dispersion is exactly f [xM^/a, It might, therefore, be suspected that '39723 is only an 
approximation to f, which may be the rigorous value. But my numerical calculations 
were carried out with so much care that I find it almost impossible to believe that 
there is an error as large as 3 in the third place of decimals, or, indeed, any error at 
all in the third figure. Moreover, it would be expected that, if this very simple 
relationship is rigorously correct, it would be possible to prove it rigorously, just as 
it is rigorously shown above that wjp = g- ; but I am unable to find any analytical 
relationships by which the minimum value of -^ yS^ can be deduced. If my arithmetical 
process be correctly carried out, then we ought to find that, when r = '6264, dyjdx^ 
should be equal to — x^d^yjdx^ or e^^jx^. Now, I find that, when r = '6264, 
dyJdxi-= Vb770?> and e^^jx^ =- l'b770, so that the two agree to four places of 
decimals. I conclude, therefore, that the true minimum of -g- yS^ is '3972.'"'' 

It will be observed that, as aja^ increases to -infinity, \ /3^ terminates by being equal 
to ^. M. EiTTEK has found that it rises above ^, and oscillates about that value an 
indefinite number of times with diminishing amplitude, gradually settling down to \ 
as a/a^ becomes infinite. The values in the preceding table are not, however, carried 
far enough to exhibit these oscillations of \ yS^. A consequence of this result is that 
there are a number of modes of equilibrium of a gas at a given temperature, provided 
that the temperature lies within certain narrow limits. This very remarkable 
conclusion is rendered more intelligible by Mr. Hill's treatment than by M. Hitter's. 

This point has, however, no bearing on. the present investigation. 

In any one of the solutions comprised in Table II. we may complete the table of 
densities by the formula (24), viz.. 



w \ey^ 



p x^^ dy^jdo(\ {x^ = a-^la) ' 



and I shall later proceed to do this in the one case which has interest for our present 
problem, namely, where the temperature is a minim^um, so that a/a^ is '6264. The full 
numerical results may be more conveniently given hereafter, and it will only be now 
necessary to indicate how they are to be computed. 

When, for example, r=: '1 X a^, rja=. •l/*6264 = '1596; thus, our equidistant values 
of the density and other functions will proceed by multiples of '1596 a up to '9578 a, 
and the limit of the isothermal sphere is where r = a. 

When the temperature is a minimum ^ j3^ = '39723, and we have Wq=: ^ p; there- 
fore, iv/wq = w/^ p, and, therefore, if yi^Q be the value of y^, when r = '626A a^, 

* This is confirmed by Mr. Hill. His equation s =^ z is equivalent to Xi^ d^y^jdx^ -f ^i dyjdx-^ = 0, and 
it appears from his tables that s = z = 2*517. Now, s = 3/y8^, and the reciprocal of 2"5]7 is '397. 

D 2 
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ivj^ p =. e^^~^^'\ Thus, for example, at the centre, ^vj^ p is 3214, and when 
T =: -4789 a it is 5'7417. 

The proportion of the mass M which is included in radius ajx is ^ /3^ dy/dx 
= ^ fi^a^ dyja dx^ = '^^i-^i- dyjdx^. Hence, the masses may be computed. 

At any part of the isothermal sphere gravity g is to be found from 

or, expressing g in terms of G gravity at the surface, we have, since G = [juM/a^, 

^ = ^/3^x^-^^, (27) 

The angular velocity of a body moving in a circular orbit at any part of the nebula, 
and its linear velocity v are also easily to be found. 



§ 5. On an Atmosphere in Convective Equilibrium, 

I shall now suppose that a sphere of gas of mass M at minimum temperature is 
bounded by an atmosphere in convective equilibrium, with continuity of temperature 
and density at the sphere of discontinuity of radius a. Let v^ be the mean square 
of velocity of agitation in the isothermal sphere, and v^ that at any other radius r. 
Then throughout the isothermal sphere v^ = v^^ but in the layer outside v^ gradually 
decreases to zero. 

Let w^ be the density and p^ the pressure at radius a, and w, p the same things at 

radius r. 

Then, if the ratio of the two specific heats be that deduced from the simple kinetic 
theory of gases, without any allow^ance for intra-molecular vibrations, we have that 
ratio equal to f . 



Hence, 



ID \3 



p^i\[z:h 



and 



also 



^0. 









Now, the equation for the hydrostatic equiUbrium of the layer is 

'^'^^--\-liM-\-iTrix^ior^dr=zO (28) 
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a v^ / tv\^ 



X = - , % 



T ' V^ \V)r ' ^ 



'0 X'^-'O/ 



and we have 



r^' dp p. •. „ dz 



47r/xa^ = - — = — , since ?^o = i p rigorously, 



^0 






Hence, our equation is 



It is obvious the -f ^^ M dzjdx is the whole mass (expressed in terms of the mass of 
the isothermal sphere) enclosed inside of radius ajx. The differential equation to be 
satisfied is 

^^ dx^^ x^"^ ^ ^ 

We have seen in the last section that ^ ^S^ = •39723, and, hence, |^^ = '99308, 
This equation is not so easy to solve as that in the last section, and I have not 

succeeded in finding the general law of the coefficients in an expansion. Nevertheless 

it is easy to find a series which will do all that is required. 

Let c be any value of x for which we know z and dzjdx, and let 

Assume 

Z '=: Zq{1 -{- A^^ + Ac^^'^ + ^3<f^ + •••}• 

Then, if the suffix indicates the value of a symbol when ^ = c and £ = 0^ we 
have 

^0 — '^O? 
dz\ __ J 



'dh^\ 

But 



^ucLo'^n* 



^dh \ 6 ^Q^ 



dxy. 5^2 c^' 
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and 

so that, if 2;q and {dzldx)^ are known, ^^ ^^ known. 

The differential equation (30) which we have to satisfy is 



5./Q2 



(Pz 



or 



1 [^ ^^'dH^^^fz 



A^\G j <ip \z, 



0; 



Now, by expansion. 



z \3 



+ ^iA + \{A,A, + \Ai)-\AM. + -hA,'W 
+ 3 \_A, + 1 {A,A, + ^,^3) - 1 {AM, + ^1^/) + -hAM, - Tk-4i^] f 

+ 3 [^e + 1 (^1^5 + ^,^, + ^Ai) - \ {A,^A, + 2A^A,A, + \Ai) 

+ rh-{A,'A, + Ui'Ai)--rhAM.+ T-h^A,^.e 
+ -ilA,^^{A^A,+A,A,+A,A,)-^A^^A, + 2A,A,A, + A,A,^+A,^A,) 

+ -h {AM, + 3^i^^2^3 + A,Ai) - xls (-i/^3 + ^AMi) 

+ 2k-4i«^2 - T(A-4 A'] 1^^ + . . (31) 

And 









5 



+ (32) 

By equating the coefficients in (81) and (32) we are able to determine the ^'s. The 
law of the series (32) is obvious, and sufficient of the series (31) is written down to 
enable us to find Ag. We can, however, obtain a good approximation to higher 
coefficients, because the coefficients in (31) become relatively unimportant. 
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^O 



We now begin the solution with 



1, % = 1, 



'dz 
dx. 







1^^ 



= 1-0070, 



'dH 
dx^ 1^ 



^ _ _ 1-0070. 



h^^^ 



Hence, 



A^ = 1-0070, -^3 = — -5035, 



whence I compute 



^3= + -41782, ^4= --30068, -^5= +'16175, ^3= --01306, ^7=~-*1333, 



^g = + -266, 



A^ — 



•378, ^10 = + '48, 



A 



11 



•6. 



With these coefficients I find 



T 



a 



z 



12 

— 11? 


12 
10 5 


12 
9 ? 


12 

8 9 


12 

7 ? 


12 

6 


— -9123 


-8160 


-7089 


•5887 


•4525 


•2982 



>. 



(33) 



Then, evaluating x~"*2J% and combining the several values by the rules for integration 
of the calculus of finite differences, I find 



|yS' 



L_ 12 


12 

10 i 


12 

9 f 


12 

8 > 


12 

7 ? 


12 

6 


dz 
dx 


1-21 


1-35 


1-527 


1-729 


1-95 



1 



V. 



J 



(34) 



When r = 2, we begin a new series with 






%n zm ^ y o iii- 



"0 



1 6^«^ 
z dx{^ 






1-9513 



•9907 X -2982 



-/xo — -— 



1 6^%' 



ZiZ CLX 










= -- 4-3686. 



+ 6-5894, 



From these I compute ^^3 = — 2*744 , ^^ = + 21*365 , A 



45-409, 



2a and 
- -094, 



^g=: + 9-932, ^7 = +-319. 

It appears that z vanishes when a; — c = — '141 or x = -359. 

It follows, therefore, that four equidistant values of x lying between r = 
r =n: a/-359 = 2*786 a correspond tox — c = 0, jr — c=-— -047, £;c — c = 

X "~" C — — *""" 1 41: 1 • 

For the first of these, where r = 2 a, we have z = -2982, and for the last, where 
r = 2-786 a, 2; = ; and, when cc — • c == — -047, or r == a/-453 = 2*208 d^ I find 
z = -2031 ; and, when x — c = — '094, or r = a/*406 == 2*463 a, I find z = -1033. 
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Finding x''^?} for these four values and combining them by the rules of integration, 
I find 

I ^ ^ = 2-1767, when r - 2786 a. ..... (35) 

We thus see that the mass of the whole system is 2*1767 times the mass of the 
isothermal nucleus, and its radius is 2*786 times the radius of the nucleus* 

The mass of the isothermal nucleus is thus 46 per cent, of the whole. M. Ritter, 
taking the ratio of the specific heats as f instead of f , says that the proportion is 
about 40 per cent. 



§ 6. On a Gaseous Sphere in ^^ Isothermal- Adiahatic^^ Equilibrium, 

M. RiTTER calls a sphere, with isothermal nucleus and a layer in convective 
equilibrium above it, a case of isothermal-adiabatic equilibrium. Since the height of 
an atmosphere in convective equilibrium depends only on the temperature at the base, 
and since the isothermal nucleus in our numerical example is at minimum temperature, 
the thickness of the adiabatic layer is a minimum, and the isothermal nucleus a 
maximum. 

We are now in a position to collect together all the numerical results of the last 
two sections in a form appropriate for our subsequent investigation. It will be 
convenient to refer all the densities and masses to the raean density and mass of the 
isothermal nucleus. Gravity may also be referred to gravity G at the limit of the 
isothermal nucleus, and velocity to v^^ the mean square of velocity of agitation in the 
isothermal nucleus. 
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§ 7. On the Kinetic Energy of Agitation and its Distribution in an Isothermal- 

Adiahatic Sphere of Gas, 

We shall now consider what would be the distribution of kinetic energy in the 
nebula if each meteorite (or molecule) were to fall from infinity to the neighbourhood 
where we find it, and were to retain that energy afterwards. This will give the dis- 
tribution of energy in a swarm of the supposed arrangement of density, if the rate of 
diffusion of kinetic energy were to be infinitely slow, and if there were no loss of energy 
through imperfect elasticity. 

The square of the velocity of a satellite in a circular orbit is one half of the square 
of the velocity acquired by the fall from infinity to the distance of the satellite from 
the centre. If the concentration has proceeded as far as radius r^ and if a meteorite 
falls from infinity to distance 7% then, if U be its velocity, and v the velocity in a 
circular orbit at distance r, 

^U^" ziz v^ =. ^ /3^ '-— . cc -r- = -J v^^x -—, in the isothermal sphere, 

Cu CLOG CLOG 






^ ^ — . a? — = f Vq^x 'y , in the adiabatic layer. 



In these formulae, by the definitions of y and z, 

-^ j in the first, and 2; = ( — j in the second. 

From these formulae v^ was computed in Table III. The value of v^ or ^U^ gives 
what may be called the theoretical value of the kinetic energy, because it gives us a 
measure of the amount of redistribution of energy by diffusion and loss of energy 
by imperfect elasticity, which must take place before the whole system can assume the 
form of an isothermal adiabatic sphere. 
. We will now go on to consider the total potential energy lost in condensation. 

We have seen that the potential energy lost by the fall of a single meteorite is 
■^Vq^x dy/dx in the isothermal part, and §Vq^x dz/dx in the outer part. 

Now, in the isothermal part a spherical element of mass is 



do(? 



~W^'-^idx, 



and the energy lost by its fall is 



^iVLp t^Q , X .J ^ g ax. 
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Hence, the whole energy lost in the concentration of the isothermal nucleus is 



But 






dx d%^ 



,00 



dy dhj 7 

iAJ ^ 7 ^-^*^ 

"I d/JU Ka/iXj 



/.CO 



ey dy 
Y^ dx 

3 i^A ■ 

\(XXJn 



f>(Xi 



c^a? = 3 



^ 



x^ 



dm — 6^« 



e^« 



9^ 

/3^ 



^w. 







/8V 



m 



Hence, the energy lost is Mv^ f 1 ~\ But in an isothermal sphere of minimum 

temperature Wq = |-p, and thus the total lost energy is ^Mv^. 
Again, in the adiabatic layer an element of mass is 

- IM^ ^^dx= + M^-, dx, 



d^x^ 



X 



and, therefore, the energy lost by its fall from infinity is 

^ ^ x^ dx 

and the whole loss of energy is the integral of this from x = 1 to i;^ = *359. When 
cc = 1, 2: = 1, and when x = 'SSO, z = 0, Hence 



rl 



z^ dz 1^^ _. 2. I 6 
— 5 t" 5 



■359 '^ ^^ 



dx 



1 ^ 



—7 (A/X» 

•359 '^ 



Thus, the whole energy lost in the adiabatic layer is 



Mv 



3 







3 1* 



•859 



5. 
~X^ 



d\ 



X 



Add this to the energy found before for the isothermal part, and the whole lost 
energy of the system is found to be 



Mv^^ 



I + 



-4 dx 

J .359 X 



* • • • • • • • lOOl 



NDwlet us evaluate the total kinetic energy existing in the form of agitation of 
molecules. 

E 2 
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In the isothermal part it is clearly \Mv^, In the adiabatic part it is half the 
element of mass into the square of velocity of agitation integrated through the layer, 

that is to say, \ . M~^dx X v^^ and, since ^ =: — , we have 



Wvn' 



+ ^, dx 



'369 ^^ 



for the total internal kinetic energy of agitation. This is rigorously one-half of the 
energy lost in concentration. 

Hence, if a meteor swarm concentrates into this arrangement of density, one half of 
the original energy is occupied in vaporising and heating parts of the meteorites on 
impact, and the other half is retained as kinetic energy of agitation. 

,^dx =^ •643. Hence, the potential energy lost in 

•359 ^ 

concentration is Mv^ {I'Qid), and that pai't of it which is retained as energy of 
agitation is ^i¥t?Q^(l'643). The whole mass of the system is 2*1767 i¥, and we may, 
therefore, write these 

•7548 (2-1767 ikf)t;o2 and | X '7548 (21767 If) V- 

It is clear then that the average mean square of velocity of agitation of the whole 
system is *7548 v^,^^ Or, shortly, the average temperature is very nearly f of the 
temperature of the isothermal nucleus. 

It follows from this whole investigation that for any given mass of matter, arranged 
in an isothermal-adiabatic sphere of given dimensions, the actual velocities of agitation 
are determinable throughout. 

§ 8. On the '^^ Sphere of ActionJ' 

When two meteorites pass near to one another, each will be deflected from its 
straight path by the attraction of the other. The question arises as to whether the 
amount of such deflection can be so great that the passage of two meteorites near to 
one another ought to be estimated as an encounter in the kinetic theory. 

We shall now, therefore, find the deflection of two meteorites, moving with the 
mean relative velocity, when they pass so close as just to graze one another. 

The mean square of relative velocity in the isothermal portion is 2t;^Q, and this may 
be taken as the square of the velocity at infinity in the relative hyperbola described 
The angle between the asymptotes of the hyperbola is the deflection due to this sort 
of encounter. 

Let ce, € be the semi- axis and eccentricity of the hyperbola. Then, if ebe large, the 

* M. RiTTEE gives '741 in place of '755, but, as already remarked, lie uses a different value for tlie 
ratio of tlie specific lieats. 
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angle between the asymptotes is 1/e ; and, if ^s be the radius of either meteorite, 
the pericentral distance (when they graze) is s. Therefore, 



5 = a (e — • 1). 



By the law of central orbits 



Therefore, 



2Vo 









But, since v^ = ^^ [jbM/a^ we have 



Ms 

e = 2;8^ — +1. 

ma 



The unity on the right-hand side is negligible, and, since 180/'n"€ is the deflection in 
degrees, that deflection is 

180 ma , 

degrees. 



27r^2 Ms 

N0W3 if 8 be the density of the body of a meteorite, m = ^ wSs^^ and, therefore, this 
expression becomes 

13^ ^ M' 
Let us find what s must be if the deflection is 10^ ; we have 



s = ^a/ 



2M 



We may, for a rough evaluation, take ^S as unity instead of \/6/5, and suppose a to 
be equal to the distance of Neptune from the Sun (viz., 4 '5 X 10^* cm.), and, as a 
very high estimate of the value of S, let us suppose the density of a meteorite is 10. 
Then, since the Sun, M^ = 2 X 10^^ grammes, and M is about a half of the Sun'g 
mass, we have 

^ "^ "^^ =(15 X 10^6)1^. 4 X lOl 



s = 



3 X 4-5 X 10^^ X 10 



Hence, m = ^Trhs^ = ^7rX 10 X 64 X 10^*= 3 X 10^^ grammes, in round numbers. 

But the Earth's mass is 6 X 10^*^ grammes, and therefore the meteorites are one- 
twentieth of the mass of the Earth. 

It follows, therefore, that, with such small masses as those with which the present 
theory deals, the deflection due to gravity is insensible, and we need only estimate 
actual impacts as encounters. 

Hence, the radius of the sphere of action of a meteorite is identical with the 
diameter of its body. 
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§ 9. On the Criterion for the Applicability of Hydrodynamics to a Swarm of Meteorites. 

The question at issue is to determine within what hmits the quasi-gas formed by a 
swarm of coUiding meteorites may be treated as a plenum, subject to the laws of 
hydrodynamics. The doctrines of the nebular hypothesis depend on the stability of 
a rotating mass of fluid, and that stability depends on the frequencies of its gravita- 
tional oscillations. Now the works of PoiisrcARE and others seem to show that 
instability, at least in a homogeneous fluid, first arises from one of the graver modes 
of oscillation, and the period of the gravest mode does not differ much from the 
period of a satellite grazing the surface of the mass of fluid. Then, in order that 
hydrodynamical treatment should be applicable for the discussion of such questions of 
stability, the mean free time between collisions must be small compared with the 
period of such a satellite. Another way of stating this is that the mean free path of 
a meteorite shall be but little curved, and that the velocity of a meteorite shall be but 
little changed by gravity in the interval between two collisions. This must be fulfilled 
not only at the limits of the swarm^ but at every point of it. The condition above 
stated will be satisfied if the space through which a meteorite falls from rest, at any 
part of the swarm, in the mean interval between collisions is small compared with the 
mean free path. If this criterion is fulfilled, then, in most respects which we are 
likely to discuss, the swarm will behave like a gas, and we must at present confine 
the consideration of the matter to this general criterion. 

It would be laborious to determine exactly the space fallen through from rest, 
because gravity varies as the meteorite falls, but a sufficiently close approximation 
may be found by taking gravity constant throughout the fall and equal to its value at 
the point from which the meteorite starts. 

We have already denoted by g the value of gravity at any part of the swarm, and 
have tabulated it in Table III. in terms of G or iiMja^. 

Now the mean interval is T= L/(v^8/3it), Hence, if D be the distance fallen in 
this time, 



-^ — 2 bf^ — 2 ...3 Q 

But 



2 8 



= 11^^'"^ and ^ 



Therefore, 

^ — _L . :Mo . !?! . r^ . bl^l 1 _ i_ . :^ . ?!?: 



'Fi' ... (37) 



The factor F^ has been tabulated above in Table III., and it increases from the 
centre to the outside. 

This criterion may be regarded from another point of view, for, if the meteorite be 
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describing a circular orbit about the centre of the swarm, D is the deflection from 
the straight path in the mean interval between two coUisions. Then the criterion is 
that the deflection shall be small compared with the mean free path. 

We may consider the criterion from again another point of view, and state that 
the arc of circular orbit described in the mean interval shall be a small fraction of 
the whole circumference. 

The linear velocity v in the circular orbit is given by 

And the mean interval T = L/\y ^^ S/Stt], Hence, if A be the arc described with 
velocity v in time T, 



^^ Zir L^ a P' a ^ . 3'/. ^^^ 

j£^ — ^--~. ~^^~-~ . -^- = ----- . '^^ nearlv, smce ;r-^ = '988. 



8^3 ' ^3/^3 • g^ ^ #/^,3 • Q^ nearly, smce ^^^ 



>7r 
: • "- nearly, smce ^ 

But the whole arc of circumference C is 27Ta/x, 

Therefore, 

A LI [glGr]^ ,^^ 



C 2a IT vjv, 







I Mq m 



{gjGfx^ 1 ____ I M^ m , 

X'o* • » » COol 



The factor F^ has been tabulated above, in Table III. 



§ 10. On the Density of Meteorites and Numerical Application, 

It is necessary to make assumptions both as to the mass and the density of the 
meteorites. We have a right to assume, I think^ that the density S is a little less 
than that of iron, say about 6, and we may put ^ttS equal to 25. Then we have 

m = IitSs^ = -8%^, and -^ = -^s. 

There is but little information about the average size of meteorites ; but, if we retain 
the symbol s, it will be easy, by merely shifting the decimal point in the final results, 
to obtain results for all sizes. Thus, if 5 = 1 cm., m = 3|- grammes ; if 5 = 10 cm., 
m = 3^ kilogrammes; if 5 = 100 cm., m = 3|- tonnes, and if s = 1000 cm., 
m= 3125 tonnes. T shall, therefore, keep s in the analytical formulae, and put it 
equal to unity in the numerical results. 

In the first place, making no assumptions as to the density or masses of the 
meteorites, we have 

M^=2'1767 K M, i^^= -39723. 
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Then, by substitution in (10) and (11), we have 



Vq = 



U X 109-86918-10 



L=lX 100-83781 



T=T X 10 



0-46863 






B 
L 

A 



I X 100-88781 X F, X *"■ 



2a 

^, = f X 100-83781 XF^X 



m 

o3 



">k 



(39) 



We will now apply this solution to a case which will put the theory to a severe 
test. Suppose that the limit of the sphere of uniformly distributed energy of 
agitation is nearly as far as the planet Uranus, so that, say a = 16aQ. Then the 
extreme limit of the swarm is at 44|aQ ; but the orbit of the planet Neptune is at 
30aQ, so that the limit is further beyond Neptune than Saturn is from the Sun. 
Now, if a/a^ = 16, I find 

u = ioS"86ii7 QYa, per sec. 
= lO'-^^'^^^aQ per annum, 
r = 10^'*6808 seconds 

= 106-96899-10 ye^^g 

Introducing these values in (39) and putting ~^s for m/s^, I find 



(40) 



V 

L 



1 '141^0 per annum = 5*374 kilom. per sec. 

s 



•9540-10 



a, 







107'9335-10 ^ 



T = 



[^/i-p] 

S 



[vlv^] [wiy] 



r 



(41) 



D 
L 

A 
C 



^ IQQ'md-io^jp^ 



=: 1 06-4489 --IO^.JT; 



Now we have in Table III. the logarithms of the several factors, which occur last in 
these formulae (41), at various distances from the centre. 

It will suffice for our purpose only to take every other value from Table III. The 
distances from the centre are expressed in terms of the astronomical unit distance, 
viz., the Earth's mean distance from the Sun. The mean free path is expressed both 
in ihe same unit and in kilometres ; and the mean intervals between collisions in 
days. The criteria D/L and A/C are, of course, pure numbers. Table IV., as it 
stands, is applicable to meteorites weighing 3-g- grammes, bat by shifting the decimal 
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point one place to the right in the last four rows of entries it becomes kilogrammes, 
one more and it becomes tonnes, and another, thousands of tonnes, and so on. 



IV. — Table of Results. 

The meteorites weigh 3|- grammes, and have the density of iron, 
extends to 4:4|aQ, a^ being Earth's distance from Sun. 



The swarm 



Distance froml r _ 
cen^tre J a^ 


Sun. 


Asteroids. 


Saturn. 


16 


Uranus. 


24 


IN'eptune. 


441 


oo 

oo 

oo 

. 00 
00 





2^55 


' 7^66 


12-77 


19-2 


32 


Velocity of Tnean") 
square in kilo- ?v =: 
metres per sec. J 


5-37 


e5-87 


537 


5-37 


5-37 


4-85 


4-12 


2^93 


Mean free path, — = 
L kilom. — 


•00028 
41,600 


•00038 
57,000 


•00157 
233,000 

•545 


•00511 

760,000 


•00900 
1,340,000 


•0122 

1,810,000 


•0199 
2,960,000 


•0552 
8,210,000 


Mean free time, 1 ^ 

in days j ~ 


•097 


•133 


1-78 


3*13 


4-70 


9^02 


35-17 


Criterion, 

L 


• • 


•0000167 


•0000832 


•000195 


•000278 


•000387 


•000709 


•00279 


Criterion, — 


• • 


•0000113 


•0000295 


•0000633 


•0000895 


•000111 


■000174 


•000497 



The incidence of the several planets in the scale of distance is roughly indicated by 
the names written above. 

The criteria shov^ that, if the meteorites v^eigh 3 §• kilogrammes, the collisions are 
frequent enough, even beyond the orbit of Neptune, to allow the kinetic theory of 
gases to be applicable for such problems as are in contemplation. For, when r/a = 32, 
the two criteria, (with decimal point shifted one place to the right) are '028 and *005, 
both small fractions. But, if the meteorites weigh 3| tonnes, the criteria cease to be 
very small, about r/a = 24. If they weigh 3125 tonnes, the applicability will cease 
somewhat beyond where the asteroids now are. 

I conclude, then, from this discussion that we are justified in applying hydro- 
dynamical treatment to a swarm of meteorites from which the solar system originated, 
even in the earhest stages of the history of the swarm. 

This discussion has, of course, no bearing on the fundamental hypothesis that 
meteorites can glance from one another on impact with a virtually high degree of 
elasticity; nor does it do anything to justify the assumption that a swarm will consist 
approximately of a quasi-isothermal nucleus with a quasi-adiabatic layer over it. 
This latter assumption I have been led to by the considerations to which we now pass. 
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§ 11. On the Diffusion of Kinetic Energy and on the Viscosity. 

In order to discuss these questions, it will be well to begin with a simple case of 
fluid motion. 

Consider two-dimensional motion, in which there are a number of streams of equal 
breadth moving parallel to y with velocity F, and, interpolated between them, let 
there be strata of quiescent fluid ; suppose then that we wish to find the motion at 
any time after this initial state. Let the boundaries of the streams V be from x = ml 
to I (2m +1)1, Then, if u be the velocity at x, parallel to y at time t, and v the 
kinetic modulus of viscosity, the equation of motion is 

clu dht 

dt dx^ 

The solution of this being of the form e"^'* cos^p^r, the complete solution satisfying the 
initial condition is — 



2V 



e-'^'^'i'' cos "^ - i e.-'-'^'i'^ cos ^ + \ e-''^'^'f'' cos ^^ 

L L V 



Now, if v/e refer time to a period r, where r = V'lrr'v^ then after a time ^r, which is 
greater than r, the solution is sensibly 



u^\V 



^ , 4 irx 

1 + — , COS -y- 



It is clear that the maximum of %i occurs when ir := 0, and the minimum when x:=l, 
and that they are 



1 V 

2 ^ 



. 4 
1 ±-~, 

ire 



Hence, the difierence between the maximum and minimum is 4 VJTre^, Therefore, 
the ratio of the greatest difference of velocities after time Or to the initial difference of 
velocities is 4/7re^ When 6 is 1, 2, 3, this ratio assumes the values 1/2'135, l/5'804, 
1/15*73 respectively. Thus, after three times the interval r, the difference of velocities 
is small. The time r may be therefore taken as a convenient measure of viscosity. 

In our problem the streams must be taken of a width comparable with the linear 
dimensions of the solar system. I therefore take Z, the width of the streams, as equal 
to (Xq, the Earth's distance from the Sun, and we have 



- » 



Now, according to the kinetic theory of gases, the kinetic modulus of viscosity is 
I/tt into the mean free path multiplied by the mean velocity. Hence, 
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V 

Hence, we have 






%\ 1% 



L \v 'V 8 



Vt- •••••• • (42) 



If we apply this formula to the solution which has been already found in Table IV., 
we obtain the following results : — 

- 0, 2-55, 7-66, 1277, 16, 19-2, 24, 32, 

T years 1082, 792, 193, 59*2, 337, 27-5, 19-8, 617. 

These results are applicable to meteorites weighing 3|- grammes in a swarm 
extending to 44| a^. If the meteorites weigh 3| kilogrammes, the values of r would 
be one-tenth of the tabulated values. If the streams were ten times as broad, the 
periods would be a hundred times as long. 

Now the periods r in the above table, even if multiplied by a thousand, must be 
considered as short in the history of a stellar system. It thus appears that the quasi- 
viscosity must be such that a swarm of meteors must, if revolving, move nearly 
without relative motion of its parts, at least in the early stages of its evolution. 

But let us consider the values of r at different epochs in the history of the same 
system. If a be the radius of the isothermal sphere the formulae (9) and (10) show 
that L/a^ varies as a^, whilst v/a^ varies as a~^. Hence r varies inversely as a\ 
Thus, as the swarm contracts, the periods r increase rapidly. 

Thus, later in the history, the viscosity will probably fall off so much that equalisation 
of angular velocity may be no longer attained, and we should then have the central 
portion rotating more rapidly than the outside, with a gradual transition from one 
angular velocity to the other. 

The modulus v gives, besides the viscosity, the rate of equalisation of the kinetic 
energy of agitation ; this corresponds in a true gas with the conduction of heat. The 
conclusion at which we thus arrive appears to justify the assumption that the whole 
of the central part of the swarm is endued with uniform kinetic energy of agitation, 
and that the mass of the quasi-isothermal nucleus is the greatest possible. With 
regard to the assumption that the nucleus is coated with a layer in adiabatic or 
convective equilibrium, it may be remarked that the velocity of agitation must 
decrease when we get to fche outskirts of the swarm, and convective equilibrium will 
probably satisfy the conditions of the case better than any other. Further considera- 
tions will be adduced on this point in the Summary. 

* Meyer, * Kinetisclie Theorie der Gase,' p. 321, The I/tt is derived from a numerical quadrature 
wliicli gives tlie value "SIS, and it is apparently only accidentally equal to I/tt. Tlie v \/ (S/Stt) is the 
mean velocity denoted Q by Meter. 
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§12. On the Rate of Loss of Kinetic Energy through Imperfect Elasticity, and on the 

Heat Generated, 

In a collision between two meteorites the loss of energy is probably proportional to 
their relative kinetic energy before impact. Therefore, the amount of heat generated 
by a single meteorite per unit time is proportional to the kinetic energy (say h) and to 



2 



m 



the frequency of collision. By (10) the frequency^ or reciprocal of T, varies as vws 
but mh is equal to (2A)% and s^mT'^ varies as m""«. Hence, the frequency of collision 
varies as h^wm"^, and the amount of heat generated by a single meteorite per unit 
time varies ^^WwnrK But, if j9 be the quasi-hydrostatic pressure, p varies as hwm"^, 
and, therefore, the heat generated by a single meteorite varies as hpw/. 

Then, to find the total heat generated per unit time and volume, we have to multiply 
this by the number of meteorites per unit volume, that is to say, by ivm"^-, which is 
equal to 3pA"*^. 

Thus the amount of heat generated per unit time and volume is proportional to 
p^m/h"^. With meteorites of uniform size, and with uniform kinetic energy of 
agitation, this becomes simply the square of the hydrostatic pressure. 

The mean temperature of the gases volatilised by collisions must depend on a 
variety of considerations, but it would seem as if the temperature would follow, more 
or less closely, the variations of heat generated per unit time and volume. 

§ 13. On the Fringe of a Sioarm of Meteorites, 

The law of distribution of meteorites found above depends on the frequency of 
collisions. But at some distance from the centre collisions must have become so rare 
that the statistical method is inapplicable. There must then be a sort of fringe to the 
swarm, which I attempt to represent by supposing that beyond a certain radius a (not 
the same as the former a) collisions never occur, and each meteorite describes an orbit 
under gravity. 

Now, at any point gravity depends on the mass of all the matter lying inside a 
sphere whose radius is equal to the distance of that point from the centre of the 
swarm. Hence, the value of gravity depends on the law of density of distribution of 
the meteorites, which is the thing which we are seeking to discover. 

We suppose, then, that from every point of a sphere of radius a a fountain of 
meteorites is shot up, at all inclinations to the vertical, and with velocities grouped 
about a mean velocity, according to the exponential law appropriate to the case. As 
many meteorites are supposed to fall back on to the surface as leave it, and this 
inward cannonade against the boundary of the sphere exactly balances the quasi- 
gaseous pressure on the inside of the sphere. Thus, the ideal surface may be annihi- 
lated. Since the falling half of the orbit of a meteorite is the facsimile of the rising 
half, we need only trace the body from projection to apocentre, and then double the 
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distribution of density which is deduced on the hypothesis that all the meteorites are 
rising. Agahi, since every element of the sphere shoots out a similar fountain, and since 
collisions are precluded by hypothesis, we need only consider the velocity along the 
radius vector. As far as concerns the distribution of density, it is the same as if each 
element shot up a vertical fountain ; but, of course, in determining the vertical 
velocity, we must pay attention to the inclination to the vertical at which the 
meteorite was shot out. 

The mass of the matter inside the sphere, v^hose attraction affords tl e principal 
part of the force under which the meteorites move, is siy M, and, for the sake of 
simplicity of notation, we shall take ^ixMja as being unit square of velocity. 

Now, let ^ (^ (r) be the potential at the point whose radius Is r, and suppose that a 
meteorite Is shot out from a point on the sphere with a velocity u, and at an Inclina- 
tion € to the vertical ; then, if r, 6 be the radius vector and longitude of the meteorite 
at the time ^, the equations of conservation of moment of momentum, and of energy 
are — 

o dd 
r^ -~ =: ua sm €, 
at 



(sT + (' ST - * (••) = «= - * («)• 



If we write f{r) =. <^ (a) — <^ (r), and eliminate dO/dt, we get 

Now, we are to regard dr/dt as the vertical velocity In a fountain squirting up from 
a point on the sphere. Then, since f{cc) =■ 0, It follows that at the foot of the 
fountain r^ dr/dt Is equal to a^u cos e. If, therefore, § be the density at the height r, 
and Sq at the foot, the equation of continuity is 

at ^ 

Therefore, 

S ahc cos e 



But now let us suppose that the meteorites are not only shot out at Inclination e, 
but at all possible Inclinations from 0^ and 90^. It is then clear that this expression 
must be multiplied by sin ed e, and integrated. Hence, if S now denotes the Integral 

density, 

*^o ahi cos e sin e d e 



8=C 



r {r^ (n^ ~~ /(^) ) ~~" '^^^^^^ ^i^^" ^V 
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where C is a constant which it will be unnecessary to determine, and where the limit 
€q will be the subject of future considerationo 
Effecting the integration, we have 



8 = 



o 



UT 



{r^ {v? — y (r) ) -— u^u^ sin^ e}% between limits, 



C 



ur 



[ {r^ {ii^ "~"/W ) Y "^ {^^ (^^ — /(^) ) — '^'^^^^ sin^ ^o)^]- 



It is obvious that, if u^ is greater than r^f(r)/{r^ — a^), the square root involved in 
dr/dt does not vanish for any value of e ; and, hence, we must simply take e^ = 90°. 
If, on the other hand, i^^ is less than this critical value, €q is that value of e which 
makes dr/df vanish. 

Thus, our formula divides into three, viz. : — 

fir) ^ 



1st. # OTeater than 



S = 



1 

C 
u 






a 

n.3 



2\l 



U 



2 



f(r) 



a^jr'^ 



2nd. u^ less than 



/« 



2/^3 ' 



cr/r 



a 



u 



i}^-fiy)f' 



3rd. u^ less than /{r) ; S = 0. 

The physical meaning of this division is as follows : If we take a station near the 
surface of the sphere, meteorites shot out at all inclinations, even horizontally, reach 
the height of our station ; and, when they are shot out horizontally, e = 90^ If we 
go, however, to a higher region, there is a certain inclination which just brings the 
meteorites at apocentre, where drjdt = 0, to our height ; but those shot out more 
nearly horizontally fail to reach us. Still higher, not even a meteorite shot up 
vertically can reach us^ and the density vanishes. 

These results only correspond to a single velocity u ; but, if v^ be the mean 
square of the velocity, the number of meteorites whose velocities range between u and 
u-^-duh proportional to u^e'"^'''^'^''' du.^ Hence, we have to multiply S by this expres- 
sion, and integrate from u ■= oo to u 2=z 0. 

Now, the first term of the first form for S is the same as the second form ; and in 
the third form 8 is zero ; hence, this first term when multipHed by the exponential has 
to be integrated from it^ = o) to/(r). The second term of the first form of 8 has to 



be multiplied by the exponential, and integrated from u^ 
Now, for the first term put 



CO toy(r)/(l — a^/r^). 



2v^ 



(^^ -/(-)) 



X 



2 . 



* OsKAR Meyer, 'Die Kinetisclie Theorie der Gase/ 1877, pp. 271-2. 
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therefore, 

u (u^ — /(^'))'' du = i^] x^ dx, 



and the limits of x are oo to 0. 
Hence, the first term is 



^Y2l'Yg-8/(r)/2.^| ^^e-^'dx 



Again, for the second term put 



3 / . f{r) 



2v^ \^^'^ 1 - a^lr^; ~ ^^^ 



and similarly introduce it into the second term, and we have 



cCyJU - ^'Ye3/(^wa-.w 



/.oo 

x'" e 





^ ^ ^ dx. 



From these expressions we may omit the constant factors ; and, if w be the density 
at height r, whilst lu^ is the density at the sphere, 






In this formula unit square of velocity is 'liiMja ; but we have elsewhere written 
v^ = ^^[jiM/a ; hence, if the special unit of velocity be given up, we may write fi^ in 
place of 2v^, and the result becomes 



9\ 1 



- = e"^A>/^^ _ 1 ^ ^ ^-3/WMi-«w ..... (43) 

It is interesting to observe the connection between this law of density and that 
which would have held if the gaseous law (due to collisions) had obtained. In that 
case, since ^(j) (r) is the potential, we should have had 

w dr ^ dr ' ^ ^ 

Now, p = ^v^w, and, therefore, 



3 



Thus, 



or 





3 

- log w^ - 2^3 (^ (a). 


w 


- 2,^ [<^ (^) - ^ (^)] ^'= -^ ^ /(^^)^ 
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The first term of our result, then, is exactly that resulting from the gaseous law, 
and the second subtractive term represents the action of the diminished velocity with 
which the meteorites move in the higher regions, when they are liberated from the 
equalising effects of continual impacts. 

By previous definition, fjiMf{7')/a is the excess of the potential at radius a above 
its value at radius r ; hence, 



fM^f 



a 



/('•) = 



4r^fJb 



J« r 



3 



wr^ dr . dr. 







Now, since /"(r) is only required for values of r greater than a, we may put iv equal 
to its mean value p, between the limits and a. Thus, 



Hence, 



rr 







wr^ dr = lur^ dr + \pa?. 



J a 



/(^) = 



Anra 



'r 1 



J a 



rpi 



wr^ dr + spci' 



a 



'^ dr 



J a 



(ViA 



= 1 



€1 



+ 3 



[rla \ 



W 



1 ^'^hP 



'?'dz. 



If this form for f{r) were substituted in (43), we should obtain a very complicated 
differential equation for w. We may, however, find two values of f{r) within which 
the truth must lie. 

First, if we neglect the attraction of all the matter lying outside of radius a, the 
second term vanishes, and we have, 



/{r) = 1 



and the law of density is 



a 
r 



w 



^-8(l-.7/r)/33__ / -|^ 



a 



0\ 1 



g8/(l+«/r)/3^^ 



(44) 







Secondly, we may suppose the density to go on diminishing according to the inverse 
square of the distance. We have seen in the preceding solution and tables that this 
is roughly the law of diminution for a long way outside the isothermal nucleus. 
According to this assumption, w == iVq a^/r^. Hence, in the second term of f[r) we 

put IV = iVq a^/r^- = Wq/z^. 
Hence, 



and 



f 



IP P h p 



1). 



/{r) = 1 



a 



+ '^r'^dz 



Z' 



ip) \ 



1 



a 
r 



+ -^ log 



T 

a 



(45) 



The substitution of this value iri (43) gives the law of density. 

In order to see the kind of results to which these formulae ]ead,]et us suppose that, when 
we have reached radius 2 in the adiabatic layer, collisions have become so rare as to be 
negligible. Then the symbols in the formulae of this section have numerical values; and, 
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in order to distinguish them, let them be accented, so that, for example, we write 
a\ yS'^, p, &C.5 in (43), (44), and (45), in place of a, yS^, p. 

Now Table III. shows that, when a = 2a, M^ == l'95M= 2M nearly. Hence, 
M'jd = Mja nearly. But, at radius 2a in Table HI., v^lv^ = •298 = *3, and this v^ 
is what we now write v^^ or fi^^fiM^/a, whilst Vq^ = fi'^[iM./a. 

But iS^ = f very nearly ; hence, fi^'^/fi^ = '3, or fi'^ = '36. 

Thus, 3/^'^ = 8-333. 

Then, substituting 2a for a, and noticing that in Table III., w/^p = *163, when 
r = 2a, the first law of density (44) becomes 



w 



iP 



= -163 



g--¥(l-3«/r)_ I I _ 






. . . (46) 



Wr 



Wf 



Again, since M' = 2M, and a = 2a, p' = ip, /, = 4 X ^^ = 4 X '1 63 by Table III., 

^ -sP sP 

and ^ = *65 = f nearly. 

Thus, according to the second assumption, we have by (45) 



3f{r) = {l 



V) + ^°^ (i) ' ^^^' ^"^'^^ ;8^2 =2-78, 






= 2-78 1 



3/« 



2a 
r 

2-78 



+ 2-78 log 



3 



2a 






2-78 



^ 



/3'2(1 - 4a2/r2) 1 + 2a jr ' 1 - 4a'V5^^ 



lo 



2(X, 



and the law of density is 



w 

Vp 



= -163 



^■r-Sf(r)l^"' 






r-^ 



. . (47) 



The values computed from these alternative formulae (46) and (47) will be compar- 
able with those in Table III, 

In Table III. we have the value of lof^p computed at distances r/a = 2*208, 2*463, 
2 '78 6. The following short table gives the result extracted from Table III. for 
comparison with the values computed from (46) and (47) : — 



r/a = 2*0, 



2-208, 2-463, 2-786. 



Convective equilib., ^ 



tp 



-163, -092, 



w 



-033, 



tP 



Second hypoth. (47), i 



w 



3P 



163, 



.A I-/ 



071, 



•029, 
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0. 



First hypoth. (46), — = -163, '074, -033, '015. 
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It appearS; therefore, that the results from the two hypotheses differ but little for 
some distance outside the region of collisions, and either line may be taken as near 
enough to the correct result. We see then that the effect of annulling collisions and 
allowing each body to describe an orbit is that the density at first falls off* more 
rapidly than if the medium were in convective equilibrium, and that further away the 
density falls off less rapidly. At more remote distances the density would be found 
to tend to vary as the inverse square of this distance. Thus, the formulae ^vould make 
the mass of the system infinite. In other words, the existence of meteorites with 
nearly parabolic and hyperbolic orbits necessitates an infinite number, if the loss to 
the system is constantly made good by the supply. 

The subject of this section is considered further, from a physical point of view, in 
the Summary at the end. 



§ 14. On the Kinetic Theory where the Meteorites are of all sizes. 

In an actual swarm of meteorites all sizes occur, for, even if this were not the 
case initially, inequality of size would soon arise through fractures. Hence, it becomes 
of interest to examine the kinetic theory on the hypothesis that the colliding bodies 
are of all possible sizes, grouped about some mean value according to some law of 
frequency. 

If there be two sets of elastic spheres in such numbers that there are respectively 
A and B in unit volume, and if -the mean squares of the velocities of the two are a^ 
and /3^ respectively, and if a and h are the radii of the spheres of the two sets, then it 
is proved that the number of collisions between them per unit time and volume is 



2AB {a + hf [f TT {a^ + ^^) ] 



1 i'~ 



We shall now change the notation, and for a and h write s-^ and s^, and for a and fi 

write Uj and %. 

Then, if 8 be the density of the spheres, their masses are IttSs/ and firSs^^. 

The condition for the permanence of condition is that the spheres of all masses 
shall have the same mean kinetic energy. Hence, we refer the mass to a mean sphere 
of radius 9, and the velocity to a square of velocity F^. 

Then 



Si^U^^ = S.^U^ = 9^ F^. 



Thus, our formula may be written 

2AB {sj^ + So) 



s 



9 \^ . Z?^^ 



^ 



„ \^l/ \'%/ _i 



(f^n*. 



* ' The Kinetic Theory of Gases,' by H. W. Watson, p. 11. 
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But now suppose that there are spheres of all possible sizes, and that in unit 
volume the number whose radius lies between s and 6^ + ds is 



4n 



<tW' 



s^e~''l^'ds. 



y_ 



TT 



Since the integral of this from oo to is n, it follows that n is the number of 
spheres of all sizes in unit volume. 

If p be the total mass in unit volume, or the density of distribution, 



p^~ 



4% 

s/lT 

4:71 

4n 
^7~ 



J ^"^ 



00 



IttSo"^ I x^e ^^ dx 





'TT 



Za-K 



TT 



If m be the mean mass, m — pjn ; but m = ^ttS?^ ; hence, 



and 



.3 



,,3_^ /,\3 



4 



V 



TT 



•1 



'3, 



4 

\/nT 



a 



3 






a 



^3/ J 



If the A spheres of radii Si are those whose radii lie between Si and s^ + (^^i, and 
the B spheres of radii s^ are those whose radii lie between s^ and s^ + <^%> 



A 



B = 



i!LW^-..v.^^3 



V IT \(T 



(T 



4:11 /So\^ 



3 \ ^ - S//0-2 _% , 



V 



/tt \ a 



e 



cr 



Hence, the formula for collisions between the ^'s and ^'s is 



64:71^ 



TT' 



{inVJ.{s, + s,) 



2 






3 /a-\3- 

+ 



'2/ J 



'' 53_ ^ - (Si^ + SgSyo-s ^1 f^ 



e 



(T a 



or, if we write x = ''^i/cr, y = sjor, it is 



647z3 



'TT' 



{' 



lTrV^)i(r^{x-\-yf{x^-\-y^f{xyfe-^-fdxdy. . . . (48) 



* If the spheres are grouped about a mean mass, instead of about a mean radius, according to a law of 
this kind, the subsequent integrals become very troublesome. Any law of the kind suffices for the 
discussion. If, however, I had foreseen the investigation of § 16, I should not have taken this law of 
frequency. 
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But 



cr^ = V^ S^ 



and 



64 



TT' 



(f^) 



TT 



?>^ # 



O 



TT 



13 



3^ 



Hence, the number of collisions per unit time and volume between spheres whose 
radii range between s^ and s^ + ds^, and others with radii between s^ and s^ + ds2, is 



32 2« 



S oi * Vs;^"^^ . (^ + yY {x^ + 2/^)' (^y)' e'""'*'"'^' (ix (iy . 



TT^^ 3 



The number of collisions of a single sphere per unit time is l/n of this, and, since 
n = p/m, we have for the collisions of a single sphere the factor — - instead of VsV, 

Then the total number of collisions of all kinds in unit time, or the reciprocal of the 
mean free time, is the double integral of this from oo to 0. 

For the purpose of carrying out the integration^ we may conveniently, as an 
algebraic artifice, change from the rectangular axes x, y to the polar coordinates r, 6, 
Thus, 



«.00 rtOO 

(^ + yT (p^ + y^)^ i^y)^ e~^'~^' dx dy 

Jo Jo 



/»00 



fi. 



r^ e ^ dr 











TT 

(sin 6 + cos Of (1 — sin 9 cos Of (sin cos Of dO, 



Now, if we put r •=■ z^, 



00 







11 



r^ e 



'""dr 



.00 







Q f -j /.CO 

4.4,4 Jo 



For the transformation of the second integral, put 



and we find 



z = cos — sin 0, 



^(sin + cos 0)^ (1 - sin cos ^)^ (sin cos 6^)-^ d^ == f "^ '| (2 - ^^)^ (1 - z^f dz 
J ^ J - 1 



Hence, the whole integral is 







(2 - z^f^ (1 «- ^4)^ &. 



.00 



4A 
32 



e 



-^V0 



n 











(2-^^)^(l -~^*)*c^^. 



and the mean frequency of collision of a single ball per unit time is 



.i. ,^ J 



JL.5 . 
4 



3\ 2" F (29)^ 



.CO 



ri 







e-'' dz (2 -- z^f (1 — ;2^)^ c?^. 
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The second of these two integrals cannot, I think^ be evaluated algebraicallj, but 
its value is easily found by quadratures. I find, then, 

f (2 - z^f (1 — z^f dz = 1-2999. 

The former of the two integrals may be evaluated as follows :— 
Let 



then, 







e ^^ dxy 



iPs^i 



fOO /.CO 
•'0 



h Jo 



^^,4(i_X3i^.2a)^^^.^^ 






ri-TT rco 



= 2 



1 I 

77 






(1 - i sin^ 0)* 

= 77*^(45°), 

where F is the complete elhptic integral with modulus sin 45°. 
Hence, 



I=^7riFK 



1 M. 



We thus have the frequency of collision given by 



Now, Legendbe's Tables give 

logi^= -2681272, 

with which value we easily find for T the mean free time, or \jT the frequency, 

7^= 5-3318 "-^'==1^^^ (49) 

1 Tiilp mjp •^ • ^ 

If IJTq be the frequency of collision when the spheres are all of the same size and 
mass 9 and m, and are agitated with mean square of velocity F^, we have, by the 
ordinary theory, 

L^.^iL Zm^,.,,,,ZSM. (50) 

1q ^ 3 mlp m/p ^ ' 

^ I owe tliis to Mr. Forsyth, and the result verifies an evaluation by quadratures which I had made. 
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It follows, therefore, that in our case collisions are more frequent than if the balls 
were all of the same size in about the proportion of 4 to 3. 

In order to find the mean free path, we require to find the mean velocity. 

If u^ be the mean square of the velocity for any size s, the proportion of all the 
spheres of that size which move with velocities lying between v and v -{- dv m 

-~~r~- 'ifie~y^ dy, 
where y^ = 3v^/2u^. 



But the number of spheres of size between s and s + ds, in unit volume, is 

4^^ 



V^TT 



cc%""*'^ dx, 



where x = s/cr. 

Hence, the mean velocity U is given by 

16 "" "" 






vx^yh'-'^'-'y'dxdy. 



Now- 



V = y/f . uy, and shi^ = s^ F^, or x^u^ 



V^ 



4 



so that 



(T. 



\/lT 



v\ 



u 



Therefore, 



But 



2 2 /2 

I x'^^ F, and v = 1^75 x'^'-y V, 



IT 



TT 






A CO 



-co 



JQ J 







x^yh"''^'^^^ dx dy. 



therefore, 



/.OO 







3/% ^' dy =: I, and 



• 00 



CO 



FT — V "^ TT" 



00 







i^^e ^"^ dx= 2 ^^^e""^* (i^ ; 

"^O 



2:^e""^' <i^. 



This integral may be evaluated as follows :- 
Let 



,/=: 



00 







x%""^* dx., 



.CO -00 

4 e/^ =r- 4 x^yh'^'^'-^y" dx dy 



-co .1^ 



r^sin^ 2(9 e'^''^'-^"^'^'^^Wdrdd 



/«00 /.^ 



JL 

4 



O-'O 

00 /•x 



z^ sm^ (l> e.-""'-'''''"''^^ dz d<f> 



= 1 



ir 



sm' 






IT 



4^^ 



1 i 

4^ 



■J sin^ ^ 



Q (1 — -J sin^ <^) 



i d<p 



r ■ /.I 



Q (1 — l^sin^ (^)' 
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Now, 



f 



^ ___d4 E 

(1 - h^ sin^ ~ ^• 



and in the present case k^ == F^ = -|. 
Hence. 

where E and F are the complete elliptic integrals with modulus sin 45^. 
In Legendre's Tables, we find 

E =1-350644, F = 1-854075, and 2E — F = -847213. 
Then, 

~^^ ^^^(2E-F)= 1-91377. 



The mean free path 



and thus 



Z = UT =1-9138 VT=^:^^ ^, 



^ ¥m (2^* . . •. (51) 



If the spheres had all been of the same size, we should have had 

^ ■" TT (29)^ v/ 2 "^ 4-44 (29)^ . l^oz; 

Hence, finally from (49) to (52), if there be a number of spherical meteorites, of 
uniform density, of all sizes with radii grouped about a mean radius according to the 
law of error, and if S be the diameter of the meteorite of mean mass m, and p be the 
density of the distribution of meteorites in space, and ^mV^ their mean kinetic 
energy of agitation, then the mean free path L, mean free time T, and mean velocity 
U are given by 

j^ !??ve _ „5_ Wp . 

2-786 S^ ^^ S^ iicdiij'. 



5*332 FaS^ ^ ^ VS^ nearly, 
U = 1-9138 V=2F nearly. 



Ip^* • » e • . • It/Ol 



y 



Also the mean free path is about -^-j-ths, and the mean free time about f of that 
which would have held if the meteorites had all been of the same size m and had had 
the same mean kinetic energy -^ m V^, 



* I owe tMs to Mr. Forsyth. 
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§15. On the Variation of Mean Frequency of Collision and Mean Free Path for the 

several sizes of balls. 

Each size of ball has its own mean frequency of collision and mean free path, and it 
is well to trace how the total means evaluated in the last section are made up. 

We have already seen in (48) that (substituting for or its value in terms of 9) the 
number of collisions per unit time and volume between balls of sizes ^ to s -{- ds and 
balls of sizes s to s' + ds is 

^ ii^vj- ^ sM« + yf («' + yj i^yf ^-^-^ dxdy, 

where x = s/cr, y = s jcr. 

But the number of balls of size 6' to 5 + ^^ in ^^^^ volume is 

4^ o 7 

—V- x^e ^ ax. 

Hence, the mean frequency of collision for a ball of size s with all others is 

4^ • '^ (f ^^T^^ 5^) {^ + y? {^^ + y^y x-'yh-f dy. 
Now, 

Therefore, if we write l/r for the frequency of collision of a ball of size s with all 
others, we have 

1 2^7r' r(29)V?^ 



Now, the mean frequency for all sizes is given by 



{x + yY {x^ + y^)" y^ e"^' dy. 



Hence, 



1 r (2<if 

-"= 5-3318 • ^-^. 

1 mjp 



==-1780-PjJ {x^rllf{x^^y^)^r^^y'dy. ..... (54) 

The integ^^al involved here cannot in general be determined algebraically; but, if 
X be very small, or very great, we can find an approximate value for it. 
If x be very small, the integral becomes 

y'^e-y" dy—^^TT, and - = 'l 18 -- 
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If X be very large, the integral becomes 



/»00 



X' 



Now, 



foo 
z^ e " ^* dz, 




CC"' 



T 



(J 



IT 



_7_ 
12 



and 



1 



00 



^ n-S^ 



z^ e 



dz = l7T^{2E-Ff 



Therefore, the integral becomes 



S(2^-^Kr 



and with the known values of E and F this gives us 



- = •282^ 



For intermediate values of s recourse must be had to quadratures for evaluating 
the integral. I have therefore determined, by a rough numerical process, sufficient 
values of the integral to render possible the drawing of a curve for the values of 
T/r for all values of s. The following table gives the results for the integral 
ir {^ + VY (^^ + y^)^ iT ^""^^ ^y? which may be denoted by K\ — 





K 


s = is 


1-71, 


5 = |S 


2-90, 


s — s 


4-94 , 


s — is 


12^97, 


s _ 2s 


28-75 . 



If these values be introduced in the formula 



-= -1780 jrf~ ,, 



we obtain 



MDCCCLXXXIX. — A. 



H 



Th 



.9= is 


•86, 


•5- |S 


•80, 


s — S 


•88, 


s = f s 


1-26, 


s— 2s 


1-81 . 
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These values are used for forming the curve, entitled ^^ frequency of coUision," in 
fig. 1 below^ and they are supplemented by the values found above for Tjr, in the 
case where sh is either very small or very large. 

The frequency becomes infinite when the balls are infinitely small, be(jause of the 
infinite velocity with which they move, and again infinite for infinitely large balls, 
because of their infinite size; But it must be remembered that there are infinitely 
few balls of these two limiting sizes. 

We have now to consider the mean free path, say X, for the several sizes. 

If ^i^ be the mean square of velocity for the size s, the mean velocity for that size is 
u -\/(8/37r), by the ordinary kinetic theory. 

From the constancy of mean kinetic energy for all sizes, we have 

so that the mean velocity for size s is 

^(?A)V(8/37r). 

But, if U be the mean velocity, and L the mean free path, and T the mean free 
time for all sizes together, we have 

U IX 



V 



1-9138 1-9138 T 



Therefore, the mean velocity for size s is 



1-9138 U T~ ^^^Msj T 



But the mean velocity for size 5 is X/r ; hence, 

1."^" ^^^^\s) T"~" 1780 • K 
__ 2-705 

When s is very small, we find X/L = 4, and, when s is very large, X/L = 1*7 (9^)^ 
Thus, for small values of s, the mean free path reaches a constant limit 4, and for 
large values it becomes infinitely small. 

The intermediate values, sufficient for drawing a curve, are given in the following 

short table :— 

X/L. 

5 — ^9 1*58., 

s = 9 '55. 

s z= |s -21. 

ii = 2^ -09. 
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These values are set out in the annexed figure in the curve marked ''free path/' 
and are supplemented by the values found above for small and large values of s. 
The constant limit 4 falls outside the figure. The horizontal portion of the curve is 
asymptotic to the 5-axis. 

Fig. 1. 




No immediate use is made of these conclusionSj but it was proper to examine this 
point in the theory. 



§ 16. On the Sorting of Meteorites according to size and its Results, 

It is a well-known result of the kinetic theory of gases, that if a number of different 
gases co-exist J each gas has the same density as though it alone existed, and was 
subject to the resultant forces of the system ; also the mean kinetic energy of agitation 
of each gas is the same. From this it follows that the elasticity of each gas is 
inversely proportional to the mass of its molecule. 

Carrying on this conclusion to meteorites, we see that the elasticity of the gas 
formed by large meteorites is less than that for small ; and, hence, there is a greater 
concentration of large meteorites towards the centre, and there will be a sorting 
according to size. The object of this section is to investigate this point. 

In §§14 and 15, the laws of a kinetic theory were investigated when the gas 
consisted of molecules of all masses, grouped, according to a law of frequency, about 
a certain mean radius, and molecules of infinite mass were considered to be admissible, 
with, of course, infinite rarity. Now, if we were to continue to use that law of 
frequency of masses in the present investigation, we should find, as an analytical 
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result, that the mean mass in the centre of the swarm becomes infinite. The 
existence of very large meteorites in sufficient numbers to give statistical constancy 
in a volume which is not a considerable fraction of the volume of the whole swarm is 
physically improbable. We shall, therefore, treat the case best by absolutely 
excluding very large masses. When such masses occur, they must not be treated 
statistically ; this is a question which I hope to consider in a future paper. Had I 
foreseen this conclusion when the investigations of the last two sections were carried 
out, a different law of frequency of mass would have been assumed. But the results 
of those sections are amply siifficient to indicate the conclusions which would have 
been reached with another law of frequency, and, therefore, it does not seem worth 
while to recompute the results by means of a fresh series of laborious quadratures. 

Any law of frequency would suffice for our purpose which excludes masses greater 
than a certain limit and rises to a maximum for a certain mean mass. For the 
present, I do not specify that law precisely, but merely assume that at some radius, 
which may conveniently be taken as that of the isothermal sphere, where r = a, the 
number of meteorites w^hose masses lie between x and ^ + 8x is f{x) Bx ; it is also 
assumed that x may range from M ^''' to zero. 

The meteorites whose masses range from cc to x -^ Sx may be deemed to constitute 
a gas. Suppose that at radius r the number of its molecules per unit volume is Sn, 
its density Stt?^ its pressure Sp, and let the same symbols, with suffix 0, denote the 
same things at radius a. Since all the partial gases are in the permanent state, they 
all have the same mean kinetic energy of agitation, equal to |/i, suppose. Throughout 
the isothermal sphere, this h is constant, and equal, say, to h^, but varies with the 
radius in the adiabatic layer over it. It follows, therefore, that the mean square of 
the velocity of the particular partial gas a? to a? + Sx is equal to h/x, and the relation 

between Sp and Siv is 

h 

Sp =z ^~ Sm 

I^et} -— )(he the excess of the gravitation potential of the tuhole swarm at radius r 
above its value at radius a. 

Then, since each partial gas behaves as though it existed by itself, the equation of 
hydrostatic equilibrium of the partial gas a? to a? + S^ is 

Bw dr clr 

The investigation must now divide into two^ according as whether we are con- 
sidering the isothermal sphere or the adiabatic layer. 

The Isothermal Sphere, 

Here we have h a constant and equal to h^, and hp varies as hw, so that 

** This M is not to be confused with M, the mass of the isothermal sphere. 
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ha , hw 






or 

Bio Q /7 

— oxx//io^ 



Bw, 







Now it is obvious that Sn/SriQ = Siv/Swq ; and, therefore. 

But, by the definition o{f{x), 

Suq z=:f{x) hx ; 
hence, 

hn = e''^^^i^'^f{x)hx (55) 

This is the law of frequency of mass xio x ']r ^x at radius r. 
Now, if m, m^ be the mean masses at radii r and a respectively, 

rM 

X e"^^''^^^'' f {x) clx 

^=j^ ; ....... (56) 

e-3xAy(cc)(ix 

and, if we put ;)( = 0, we obtain m^ from the same formula, 

It is also clear that, if w be the total density of the swarm at radius r, 



= \xdn= \ xe'~'^^'''^''f{x)dx. ...... (57) 



By the definition of x? ^ind in consequence of the supposed spherical arrangement of 
matter, we have 

X = ~i 4:7rii%V7''^dr J dr. 

If this value were substituted in (57), we should obtain a very complicated 
differential equation to determine %Vy the solution of which is hopelessly difiicult. 
We may, however, assume without much error that the li) in the integral expressing 
y^ is the density of meteorites, all of which are of the same size m ^ and which are 
agitated with mean kinetic energy -^/^q. If this density be written w, we then clearly 
have 

The values of w and Wq may be extracted from Table III. of solutions in § 6. 
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Then we have 



3va? 



^.- 



X 



TPV 



h 







TYl 



7 log -" == qx^ suppose. 







where q is rigorously equal to — Zyjli^\ \y\xt for computing the approximate value 
(Ijm') log (w/wq) is to be employed. 

In order to proceed to the evaluation of the mean mass at various distances, we 
must assume some form iox f{x), 

T assume, then, that 



%%, 



/(.r) = ^«a;(M-^). 



M 



It is easy to show that 



fM 



dx 



n 



0? 



1 f^^ 

and - ^fi^) ^^ = iM. 
% Jo 



Hence, the mean mass m^ = ^M, and the maximum frequency is for masses equal 



to Mq, 



Then, by (56), we have for the mean mass at radius 



r 



But 







M 







x^ (M — 0?) e'^'^'dx 



m 



M 



^(M -- x)e^^dx 



1 



x^, (M — x) e^'^dx = - [6*'^ {Wq^ — 4Mg + 6) — 2 (Mg + 3)] 



2 



rM 1 

^ (M - 0?) e^-^(i^ = - [e^^^ (Mg -- 2) + (Mg + 2)]. 
Jo 9. 



/" 



• (58) 



It may be remarked that, if Mg be treated as small, we have the first of these 
integrals equal to x^-M^ (1 + f Mg), and the second equal to ^M? (1 -j- i^g), and the 
ratio of the first to the second is ^M ( 1 + "reM^?)* 

In order to evaluate m, we proceed to introduce the approximate value for q. 

Now, 

<l = -r^^SZ^ ^nd e^%:=:(-j ; 



— loff — , and e^^'^ 
m ^ Wq 



then, writing for brevity. 



P 



log f — ' 



Wq/ . ' 



we have 



w \ M/f#/ 



m 



w, 



(P^ - 4P + 6) - 2 (P + 3) 



o> 



rJf 



p* l^\^bri 



/^\ (P^2) + (P+ 2) 



. (59) 
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Also, if P be small, the approximate result is 



m 



1 j„ _i p 

J- "T" 1 0^ • 



Before proceeding to give niimerieal values for the fall of mean mass as we proceed 
outwards from the centre of the isothermal sphere, we must consider 



The Adiahatic Layer. 

In this case we assume, as before, that the ratio of the two specific heats is If, and 
we therefore have for the relatiouBhip between Sp and Zw at radius ?% 









0/ 



Hence, 



Sw .dr- 



5. 
6 






J 



cl fSiv\i 



df \ hj), 



0/ 



But, since Sj>q, Si% apply to the radius a where h == h^, a constant, 



3 %?Q h 



'0 



J 

'X 



Thus, in the adiahatic layer the equation of hydrostatic equilibrium is 



^ X dr VSwa) dr ' 



'0/ 



whence. 



X = t 



sKi^ 



X 






• • « « 



. . (60) 



0/ 



or 



8w = Sio 







1 — 



G^w 



5h 



oj 



The investigation now follows a line parallel to that taken before. 
We have Sn/SriQ = Siv/Swqj and Suq ■=if{x) dx, so that 



Sti = f 1 "- ?f Y/(cr.) hx. 



oh 







This is the law of frequency of masses lying between x and a? + Sx at radius r. 

As hn can never be negative, we see that there can be no mass greater than | Vx ' 
and, if M be the greatest positive value of the expression /(^r), there can be no mass 
greater than the smaller of f hJx or M. 
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Thus, if m be the mean mass at radius r, 






TO-^^^^ ^ 



■a 



\-^ ~ "^;)/(*) ^^^ 



(61) 



where a is the smaller of f ^q/x and M. 

If we put X = in (61), we obtain m^, the mean mass at radius a. 
It is clear also that, if iv be the total density of the swarm at radius r, 



tv = 



X dn =: 



\ ^K 



xil — -^)f{x)dx. ...... (62) 



By definition of ^, and in consequence of the supposed spherical arrangement of 
matter, we have 

^ = - f iTTfiivr'^ dr ] dr. 

If this value were substituted in (62), we might obtain a complicated differential 
equation for w. It is clear, however, that an adequate approximation may be obtained 
by assuming that the iv in the integral expressing x is the density of meteorites, all 
of which are of the same size m\ arranged in a layer in convective equilibrium, and 
with kinetic energy of agitation at the limit r = a equal to | Hq. 

If this density be written w, and if v^ be the mean square of velocity of 
agitation at radius r, we have, by (60), and in consequence of the relationship 
(w/wo)* = (v/vo)^ 



x = ibh-(^)')=wh-'' 



and 



Let 



m \ \Wq/ / " - ^ y^ 






1 _ J^ / J^ 

/3 m' \ Vq^' 



3 / ' 



for brevity ; then, adopting the law of frequency f{x) = -^x{M -- x), as before, we 
have for the mean mass at radius r 



M^ 



f 



«. / j^\f 



m = p A _„^^Z„__ 

1 X (M --' x){l — '—\ dx 

where a is equal to the smaller of M and E 



. (63) 
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The solution now becomes different according as M or /3 is the smaller. 
First, suppose M is the smaller. Then the limits of integration are M and 0. 
If we put 
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so that the numerator and denominator of m are easily integrable. 
If now we write 
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Then, since /3 = (1 — Q)/M, we have 






(64) 



This expression has a high order of indeterminateiiess when Q 
when Q is nearly equal to unity 



1, but I find that 
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Thus, the mean mass is ^M where r = a, which we know to be correct. 
Secondly, suppose that fi is smaller than M. Then effecting the integrations in the 
same manner as before, we have 
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In order to compute from the formulae, (59), (64), {66), it is necessary to make an 
assumption as to the vakie of m' the mass of the meteorites of uniform size whose 
arrangement of density is supposed to be the same as that of the heterogeneous 
meteorites. 

We have supposed that the law of frequency of masses is known at radius a, and 
that the mean mass is tbere equal to ^M. Now, inside of that radius the larger 
masses are more frequent, and outside of it the smaller masses. I suppose, then, that 
throughout the isothermal sphere m' lies half way between m^ or ^M and the 
maximum mass M, and in the adiabatic layer that it lies half way between m^ or ^M 
aud the minimum mass 0. 

Thus, inside I take m = |M, and outside m = -^M. 

As we only want to consider the general nature of the sorting process, these 
assumptions will suffice. It may also be remarked that a large variation of m is 
required to make any considerable difference in the numerical results. 

We now have — 
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Tims, our formulae are :~ 

In the isothermal sphere, from (59), 
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in the adiabatic layer. 
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The values of w/^p and of vYvq^ are tabulated in Table III., and from these 
I compute — 
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These values (together with two others in the isothermal part) are set out in fig. 2, 
and show the law of diminution of mean mass from centre to outside. 
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Diagram showing diminution of mean mass from centre to outside. 
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The evaluation of mean mass in the fringe (see § 1 3), where collisions are supposed 
to be non-existent, is not very difficult, although it involves some troublesome algebra. 
I do not give the investigation, merely remarking that it leads to almost exactly the 
same kind of law of diminution of mean mass as we have found in the adiabatic 
layer. 

§ 17. Summary. 

The first and second sections only involved arguments of a general character in 
which mathematical analysis was unnecessary. The reader who does not wish to 
concern himself with details may therefore be supposed to have passed from §§ 1 and 2 
to this Summary. 

In order to submit the theory to an adequate test, it is necessary to discuss some 
definite case of the aggregation of a swarm of meteorites, and it is obvious that the 
only system of which we possess any knowledge is oar own. It is accordingly 
supposed that a number of meteorites have fallen together from a condition of wide 
dispersion, and have ultimately coalesced so as to leave the Sun and planets as their 
progeny. The object of this paper is to consider the mechanical condition of the 
system after the cessation of any considerable supply of meteorites from outside, and 
before the coalescence of the swarm into a star w4th attendant planets. 

For the sake of simplicity, the meteorites are considered to be spherical, and are 
treated, at least in the first instance, as being of uniform size. 

It is assumed provisionally that the kinetic theory of gases may be applied for the 
determination of the distributien of the meteorites in space. -Mo account being taken of 
the rotation of the system, the meteorites will be arranged in concentric spherical layers 
of equal density of distribution, and the quasi-gas, whose molecules are meteorites, 
being compressible, the density will be greater towards the centre of the swarm. 

The elasticity of a gas depends on the kinetic energy of agitation of its molecules ; 
and, therefore, in order to determine the law of density in the swarm, we must 
know the distribution of kinetic energy of agitation. It is assumed that, when the 
swarm comes under our notice, uniformity of distribution of energy has been attained 
throughout a central sphere, which is surrounded by a layer of meteorites with that 
distribution of kinetic energy which in a gas corresponds to convective equilibrium. 
In other words, we have a quasi-isothermal sphere surrounded by what may be called 
an atmosphere in convective equilibrium, and with continuity of density and velocity 
of agitation at the sphere of separation. Since in a gas in convective equilibrium 
the law connecting pressure and density is that which holds when the gas is 
contained in a vessel impermeable to beat, such an arrangement of gas has been called 
by M. Hitter '' an isothermal-adiabatic sphere," and the same term is adopted here 
as applicable to a swarm of meteorites. The justifiability of these assumptions will 
be considered later. 

The first problem which presents itself, then, is the equilibrium of an isothermal 
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sphere of gas under its own gravitation. The law of density is determined in § 4 ; 
but it will here suffice to remark that, if a given mass be enclosed in an envelope of 
given radius, there is a minimum temperature (or energy of agitation) at w^hich 
isothermal equilibrium is possible. The minimum energy of agitation is found 
to be such that the mean square of velocity of the meteorites is almost exactly f 
(viz. ri917) of the square of the velocity of a satellite grazing the surface of the 
sphere in a circular orbit. 

As indicated above, it is supposed that in the meteor-swarm the rigid envelope 
bounding the isothermal sphere is replaced by a layer of meteorites in convective 
equilibrium. The law of density in the adiabatic layer is determined in § 5, and it 
appears that, when the isothermal sphere has minimum temperature, the mass of the 
adiabatic atmosphere is a minimum relatively to that of the isothermal sphere. 
Numerical calculation shows, in fact, that the isothermal sphere cannot amount in 
mass to more than 46 per cent, of the mass of the whole isothermal-adiabatic sphere, 
and that the limit of the adiabatic atmosphere is at a distance equal to 2*786 times 
the radius of the isothermal sphere."^ A table of various quantities in such a system, 
at various distances from the centre, is given in Table III., § 6. 

It is next proved, in § 7, that the total energy, existing in the form of energy of 
agitation in an isothermal-adiabatic sphere, is exactly one-half of the potential energy 
lost in the concentration of the matter from a condition of infinite dispersion. This 
result is brought about by a continual transfer of energy from a molar to a molecular 
form, for a portion of the kinetic energy of a meteorite is constantly being transferred 
into the form of thermal energy in the volatilised gases generated on collision. The 
thermal energy is then lost by radiation. 

It is impossible as yet to sum up all the considerations which go to justify the 
assumption of the isothermal-adiabatic arrangement ; but it is clear that uniformity 
of kinetic energy of agitation in the isothermal sphere must be principally brought 
about by a process of diffusion. It is, therefore, interesting to consider wha-t amount 
of inequality in the kinetic energy would have to be smoothed away. 

The arrangement of density in the isothermal-adiabatic sphere being given, it is 
easy to compute what the kinetic energy would be at any part of the swarm, if each 
meteorite fell from infinity to the neighbourhood where we find it, and there retained 
all the velocity due to such fall. The variation of the square of this velocity gives an 
indication of the amount of inequality of kinetic energy which has to be degraded by 
conversion into heat and redistributed by diffusion in the attainment of uniformity. 
This may be called " the theoretical value of the kinetic energy '' ; it is tabulated in 
Table III., on the line called ^'square of velocity of satellite." It rises from zero at 
the centre of the sphere to a maximum, which is attained nearly half way through 
the adiabatic layer, and then falls again. If the radius of the isothermal sphere be 
unity, then from ^ to 2 the variations of this theoretical value of the kinetic energy 

^ These results had been previously discovered by M. Rjtter. 
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are small. Since this ^theoretical value of the kinetic energy" is zero at the centre, 
there must have been diffusion of energy from without inwards, and considerations 
of the same kind show that when a planet consolidates there must be a cooHng of 
the middle strata both outwards and inwards. 

We must now consider the nature of the criterion which determines whether the 
hydrodynamical treatment of a swarm of meteorites is permissible. 

The hydrodynamical treatment of an ideal plenum of gas leads to the same result 
as the kinetic theory with regard to any phenomenon involving purely a mass, when 
that mass is a large multiple of the mass of a molecule ; to any phenomenon involving 
purely a length, when the cube of that length contains a large number of molecules ; 
and to any phenomenon involving purely a time, when that time is a large multiple 
of the mean interval between collisions. Again, any velocity to be justly deduced 
from hydrodynamical principles must be expressible as the edge of a cube containing 
many molecules passed over in a time containing many collisions of a single molecule ; 
and a similar statement must hold of any other function of mass, length, and time. 

Beyond these limits, we must go back to the kinetic theoiy itself, and in using it 
care must be taken that enough molecules are considered at once to impart statistical 
constancy to their properties. 

There are limits, then, to the hydrodynamical treatment of gases, and the like must 
hold of the parallel treatment of meteorites. 

The principal question involved in the nebular hypothesis seems to be the stability 
of a rotating mass of gas ; but, unfortunately, this has remained up to now an 
untouched field of mathematical research. We can only judge of probable results 
from the investigations which have been made concerning the stability of a rotating 
mass of liquid. Now, it appears that the instability of a rotating mass of liquid first 
enters through the graver modes of gravitational oscillation. In the case of a 
rotating spheroid of revolution the gravest mode of oscillation is an elliptic deforma- 
tion, and its period does not differ much from that of a satellite which revolves round 
the spheroid so as to graze its surface. Hence, assuming for the moment that a 
kinetic theory of liquids had been formulated, we should not be justified in applying 
the hydrodynamical method to this discussion of stability unless the periodic time of 
such a satellite were a large nuiltiple of the analogue of the mean free time of a 
molecule of liquid.* 

Carrying, then, this conclusion on to the kinetic theory of meteorites, it seems 
probable that hydrodynamical treatment must be inapplicable for the discussion of 
such a theory as the meteoric-nebular hypothesis, unless a similar relation holds good. 

These considerations, although of a very general character, will afford a criterion of 
the applicability of hydrodynamics to the discussion of the mechanical conditions of a* 
swarm of meteorites in the kind of problem suggested by the nebular hypothesis. 

* If the molecules of liquid describe orbits about one another, the analogue would probably be the 
mean periodic time of one molecule about another. 



A SWARM OF METEORITES, A¥D ON' THEORIES OF COSMOGONY. 63 

In § 9 two criteria, suggested by this line of thought, are found. They measure^ 
roughly speaking, the degree of curvature of the average path pursued by a meteorite 
between two collisions. These two criteria, denoted D/L and AjC, will afford a 
measure of the applicability of hydrodynamics in the sense above indicated. 

After these preliminary investigations, we have to consider what kind of meeting 
of two meteorites will amount to an '^ encounter " within the meaning of the kinetic 
theory. Is it possible, in fact, that two meteorites can considerably bend their paths 
under the influence of gravitation when they pass near one another ? This question 
is answered in § 8, where a formula is found for the deflection of the path of each of a 
pair of meteorites, when, moving with their mean relative velocity, they graze past 
one another without striking. It appears from the formula that, unless they have the 
dimensions of small planets, the mutual gravitational influence is practically insensible. 
Hence, nothing short of absolute impact is to be considered an encounter in the 
kinetic theory ; and what is called the radius of *' the sphere of action '' is simply the 
distance between the centres of a pair when they graze, and is, therefore, the sum of 
their radii, or, if of uniform si^ze, the diameter of one of them. 

The next point to consider is the mass and size which must be attributed to the 
meteorites. 

The few samples which have been found on the earth prove that no great error can 
be committed if the average density of a meteorite be taken as a little less than that 
of iron, and I accordingly suppose their density to be six times that of water. 

Undoubtedly, in a swarm of meteorites all sizes exist (a supposition considered 
hereafter) ; for, even if originally of one uniform size, they would, by subsequent 
fracture, be rendered diverse. But in the first consideration of the problem they have 
been treated as of uniform size, and, as actual average sizes are nearly unknown, 
results are given in the numerical table for meteorites weighing ^\ grammes. By 
merely shifting the decimal point one, two, or three places to the right the results 
become applicable to meteorites weighing 3|- kilogrammes, 3| tonnes, 3125 tonnes, 
and so on. 

It is known that meteorites are actually of irregular shapes, but certainly no 
material error can be incurred when we treat them as being spheres. 

The object of all these investigations is to apply the formulae to a concrete example. 
The mass of the system is therefore taken as equal to that of the Sun, and the limit 
of the swarm at any arbitrary distance from the present Sun's centre. The theory is, 
of course, most severely tested the wider the dispersion of the swarm; and, accordingly, 
in the numerical example the outside limit of the Solar swarm is taken at 44| times 
the Earth's distance from the Sun, or further beyond the planet Neptune than Saturn 
is from the Sun. This assumption makes the limit of the isothermal sphere at 
distance 16, about half way between Saturn and Uranus. 

The results, applicable to meteorites of 3|- grammes, are exhibited in Table IV., § 10. 

The velocity of mean square in the isothermal sphere is \/(6/5) of the linear velocity 
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of a planet at distance 16, revolving about a central body with a mass equal to 
46 per cent, of that of the Sun, viz., 5^ kilometres per second; and in the adiabatic 
layer it diminishes down to zero at distance 44^. This velocity is independent of the 
size of the meteorites. 

The mean free path between collisions ranges from 42,000 kilometres at the centre, 
to 1,300,000 kilometres at radius 16, and to infinity at radius 44^. The mean 
interval between collisions ranges from a tenth of a day at the centre, to three days 
at radius 16, and to infinity at radius 44^. The criterion D/L ranges from eo~ooo ^^ 
the distance of the asteroids, to -3 6"o~o ^^ radius 16, and to infinity at radius 44^. 
The criterion A/C is somewhat smaller than D/L, All these quantities are ten times 
as great for meteorites of S^ kilogrammes, and a hundred times as great for meteorites 
of 3|- tonnes. 

From a consideration of the table it appears that, with meteorites of 3| kilogrammes, 
the collisions are sufficiently frequent, even beyond the orbit of Neptune, to allow the 
kinetic theory to be applicable, in the sense explained. But, if the meteorites w^eigh 
3^ tonnes, the criteria cease to be very small about distance 24 ; and, if they weigh 
3125 tonnes, they cease to be very small at about the orbit of Jupiter. 

It may be concluded, then, that, as far as frequency of collision is concerned, the 
hydrodynamical treatment of a swarm of meteorites is justifiable. 

Although these numerical results are necessarily affected by the conjectural values 
of the mass and density of the meteorites, yet it was impossible to arrive at any 
conclusion whatever as to the validity of the theory without numerical values, and 
such a discussion as the above was therefore necessary. If the particular values used 
are not such as to commend themselves to the judgment of the reader, it is easy to 
substitute others in the formulae, and so submit the theory to another test. 

I now pass on to consider some results of this view of a swarm of meteorites, and 
to consider the justifiability of the assumption of an isothermal-adiabatic arrangement 
of density. 

With regard to the uniformity of distribution of kinetic energy in the isothermal 
sphere, it is important to ask whether or not sufficient time can have elapsed in the 
history of the system to allow of the equalisation by diffusion. 

In § 11 the rate of diffusion of the kinetic energy of agitation is considered, and it 
is shown that, in the case of our numerical example, primitive inequalities of distribu- 
tion would, in a few thousand years, be sensibly equalised over a distance some ten 
times as great as our distance from the Sun. This result, then, goes to show that we 
are justified in assuming an isothermal sphere as the centre of the swarm. As, how- 
ever, the swarm contracts, the rate of diffusion diminishes as the inverse f power of 
its linear dimensions, whilst the rate of generation of inequalities of distribution of 
kinetic energy, through the imperfect elasticity of the meteorites, increases. Hence, 
in a late stage of the swarm inequalities of kinetic energy would be set up ; thus, 
there would be a tendency to the production of convective currents, and the whole 
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swarm would probably settle down to the condition of convective equilibrium 
throughout. 

It may be conjectured, then, that the best hypothesis in the early stages of the 
swarm is the isothermal adiabatic arrangement, and later an adiabatic sphere. It has 
not seemed to me worth while to discuss the latter hypothesis in detail at present. 

The investigation of § 11 also gives the coefficient of viscosity of the quasi-gas, and 
shows that it is so great that the meteor-swarm must, if rotating, revolve nearly 
without relative motion of its parts, other than the motion of agitation. But, as the 
viscosity diminishes when the swarm contracts, this would probably not be true in .the 
later stages of its history, and the central portion would probably rotate more 
rapidly than the outside. It forms, however, no part of the scope of this paper to 
consider the rotation of the system. 

In I 1 2 the rate of loss of kinetic energy through imperfect elasticity is considered, 
and it appears that the rate estimated per unit time and volume must vary directly 
as the square of the quasi-pressure and inversely as the mean velocity of agitation. 
Since the kinetic energy lost is taken up in volatilising solid matter, it follows that 
the heat generated must follow the same law. The mean temperature of the gases 
generated in any part of the swarm depends on a great variety of circumstances, but 
it seems probable that its variation would be according to some law of the same kind. 
Thus, if the spectroscope enables us to form an idea of the temperature in various 
parts of a nebula, we shall at the same time obtain some idea of the distribution of 
density. 

It has been assumed that the outer portion of the swarm is in convective 
equilibrium, and therefore there is a definite limit beyond which it cannot extend. 
Now, a medium can only be said to be in convective equilibrium when it obeys the 
laws of gases, and the applicability of those laws depends on the frequency of collisions. 
But at the boundary of the adiabatic layer the velocity of agitation vanishes, and 
collisions become infinitely rare. These two propositions are mutually destructive of 
one another, and it is impossible to push the conception of convective equilibrium to 
its logical conclusion. There must, in fact, be some degree of rarity of density, and of 
collisions, at which the statistical treatment of the medium breaks down. 

I have sought to obtain some representation of the state of things by supposing 
that collisions never occur beyond a certain distance from the centre of the swarm. 
Then, from every point of the surface of the sphere, which limits the regions of 
collisions, a fountain of meteorites is shot out, in all azimuths and inclinations to the 
vertical, and with velocities grouped about a mean according to the law of error. 
These meteorites ascend to various heights without collision, and, in falling back on to 
the limiting sphere, cannonade its surface, so as to counterbalance the hydrostatic 
pressure at the limiting sphere. 

The distribution of meteorites, thus shot out, is investigated in § 13, and it is 
found that near the limiting sphere the decrease in density is somewhat more rapid 
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than the decrease correspond] Dg to convective equilibrium. But at more remote 
distances the decrease is less rapid, and the density ultimately tends to vary inversely 
as the square of the distance from the centre. 

It is clear, then, that, according to this hypothesis, the mass of the system is infinite 
in a mathematical sense, for the existence of meteorites with nearly parabolic and 
hyperbolic orbits necessitates an infinite number, if the loss of the system shall be 
made good by the supply.^ 

But, if we consider the subject from a physical point of view, this conclusion 
appears unobjectionable. The ejection of molecules with exceptionally high velocities 
from the surface of a liquid is called evaporation, and the absorption of others is called 
condensation. The general history of a swarm, as stated in § 2, may then be put in 
different words, for we may say that at first a swarm gains by condensation, that 
condensation and evaporation balance, and, finally, that evaporation gains the day. 

If the hypothesis of convective equilibrium be pushed to its logical conclusion, we 
reach a definite limit to the swarm, whereas, if collisions be entirely annulled, the 
density goes on decreasing inversely as the square of the distance. The truth must 
clearly lie between these two hypotheses. It is thus certain that even the very small 
amount of evaporation shown by the formulae derived from the hypothesis of no 
collisions must be in excess of the truth ; and it may be that there are enough waifs 
and strays in space, ejected from other systems, to make up for the loss. Whether or 
not the compensation is perfect, a swarm of meteorites would pursue its evolution 
without being sensibly affected by a slow evaporation. 

Up to this point the meteorites have been considered as of uniform size, but it is 
well to examine the more truthful hypothesis, that they are of all sizes, grouped about 
a mean according to a law of error. 

It appears, from the investigation in § 14, that the larger stones move slower, the 
smaller ones faster ; and the law is that the mean kinetic energy is the Bvme for all 

It is proved that the mean path between collisions is shorter in the proportion of 
7 to 11, and the mean frequency greater in the proportion of 4 to 3, than if the 
meteorites were of uniform mass, equal to their mean. Hence, the previous numerical 
results for uniform size are applicable to non-uniform meteorites of mean mass about 
a third greater than the uniform mass; for example, the results for uniform 
meteorites of S^ tonnes apply to non-uniform ones of mean mass, a little over 4 

tonnes. 

The means here spoken of refer to all sizes grouped together, but there are a separate 
mean free path and a mean frequency appropriate to each size. These are investigated 

[^ It must also be borne in mind that the very liigb velocities, which occur occasionally in a medium 
with perfectly elastic molecules, must happen with great rarity amongst meteorites. An impact of such 
violence that it ottght to generate a hyperbolic velocity will probably inerely cause fracture. — Added 
JSTov. 23, 1888.] 
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in § 15, and their various values are illustrated in fig. 1. The horizontal scale in that 
figure gives the ratio of the radius of each size to the radius of the meteorite of mean 
mass. The vertical scales are the ratio of the mean free path of any size to that of 
all sizes together, and the ratio of the mean frequency for any size to that of all sizes 
together. The figure shows that collisions become infinitely frequent for the infinitely 
small ones, because of their infinite velocity ; and again infinitely frequent for the 
infinitely large ones, beca,use of their infinite size. There is a minimum frequency of 
colHsion for a certain size, a little less in radius than the mean, and considerably less 
in mass than the mean mass. 

For infinitely small meteorites, the mean free path reaches a finite limit, equal to 
about four times the grand mean free path ; but this could not be shown in the figure 
without a considerable extension of it upwards. For infinitely large ones, the mean 
free path becomes infinitely short. It must be borne in mind that there are infinitely 
few of the infinitely large and small meteorites. 

Variety of size does not, then, so far, materially affect the results. 

But a difference arises when we come to consider the different parts of the swarm. 
The larger meteorites, moving with smaller velocities, form a quasi-gas of less elasticity 
than do the smaller ones. Hence, the larger meteorites are more condensed towards 
the centre than are the smaller ones, or the large ones have a tendency to sink down, 
whilst the small ones have a tendency to rise. Accordingly, the various kinds are to 
some extent sorted according to size. 

In § 16, an investigation is made of the mean mass of the meteorites at various 
distances from the centre, both inside and outside of the isothermal sphere, and fig. 2 
is drawn to illustrate the law of diminution of mean mass. 

It is also clear that the loss of the system through evaporation must fall more 
heavily on the small meteorites than on the large ones. 

After the foregoing summary, it will be well to briefly recapitulate the principal 
conclusions which seem to be legitimately deducible from the whole investigation ; 
and, in this recapitulation, qualifications must necessarily be omitted, or stated with 
great brevity. 

When two meteorites are in collision, they are virtually highly elastic, although 
ordinary elasticity must be nearly inoperative. 

A swarm of meteorites is analogous with a gas, and the laws governing gases may 
be applied to the discussion of its mechanical properties. This is true of the swarm 
from which the Solar system was formed, when it extended beyond the orbit of the 
planet Neptune. 

When the swarm was very widely dispersed, the arrangement of density and of 
velocity of agitation of the meteorites was that of an isothermal -adiabatic sphere. 
Later in its history, when the swarm had contracted, it was probably throughout in 
convective equilibrium. 

The actual mean velocity of the meteorites is determinable in a swarm of given 
mass, when expanded to a given extent. 

R. 2 
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The total energy of agitation in an isothermal-adiabatic sphere is half the potential 
energy lost in the concentration from a condition of infinite dispersion. 

The half of the potential energy lost, which does not reappear as kinetic energy of 
agitation, is expended in volatilising solid matter and heating the gases produced 
on the impact of meteorites. The heat so generated is gradually lost by radiation. 

The amount of heat generated per unit time and volume varies as the square of 
the quasi-hydrostatic pressure, and inversely as the mean velocity of agitation. The 
temperature of the gases volatilised probably varies by some law of the same nature. 

The path of the meteorites is approximately straight, except when abruptly 
deflected by a collision w^ith another. This ceases to be true at the outskirts of 
the swarm, where the collisions have become rare. The meteorites here describe 
orbits, under gravity, which are approximately elhptic, parabolic, and hyperbolic. 

In this fringe to the swarm the distribution of density ceases to be that of a gas 
under gravity, and, as we recede from the centre, the density at first decreases more 
rapidly, and afterwards less rapidly, than if the medium were a gas. 

Throughout all stages of the history of a swarm there is a sort of evaporation, by 
which the swarm very slowly loses in mass, but this loss is more or less counter- 
balanced by condensation. In the early stages, the gain by condensation outbalances 
the loss by evaporation ; they then equilibrate ; andj, finally, the evaporation may be 
greater than the condensation. 

Throughout the swarm the meteorites are partially sorted, according to size. As 
we recede from the centre, the number of small ones preponderates more and more 
and, thus, the mean mass continually diminishes with increasing distance. The loss 
to the system by evaporation falls piincipally on the smaller meteorites. 

A meteor-swarm is subject to gaseous viscosity, which is greater the more widely 
diffused is the swarm. In consequence of this, a widely extended swarm, if in 
rotation, will revolve like a rigid body, without relative movement of its parts. 
Later in its history, the viscosity will, probably, not suffice to secure uniformity of 
rotation, and the central portion will revolve more rapidly than the outside. 

[The kinetic theory of meteorites may be held to present a fair approximation to 
the truth in the earlier stages of the evolution of the system. But ultimately the 
majority of the meteors must have been absorbed by the central Sun and its attendant 
planets, and amongst the meteors which remain free the relative motion of agitation 
must have been largely diminished. These free meteorites — the dust and refuse of 
the system — probably move in clouds, but with so little remaining motion of agitation 
that (except, perhaps, near the perihelion of very eccentric orbits) it would scarcely 
be permissible to treat the cloud as in any respects possessing the mechanical 
properties of a gas.''"] 

The value of this whole investigation will appear very different to different minds* 
To some it will stand condemned, as altogether too speculative ; others may think that 

^^ Added Nov. 23, 1888. 
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it is better to risk error on the chance of winning truth. To me, at least, it appears 
that the Hne of thought flows in a true channel ; that it may help to give a meaning 
to the observations of the spectroscopist ; and that many interesting problems, here 
barely alluded to^ may, perhaps, be solved with sufficient completeness to throw light 
on the evolution of nebulae and of planetary systems. 



